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APRESENTAÇÃO 


É com grande satisfação que divulgamos mais uma organização de livro na 
temática das geociências, com o título “Descobertas e desafios nas geociências”, 
que apresenta uma série de quatro artigos com diferentes abordagens e estudos 
voltados para estudos de análise ambiental em diferentes contextos. 

Além da abordagem ambiental em áreas de bacia hidrográfica no estado 
de Minas Gerais até estudos para análise de modelagem geoquímica aplicada. 

Assim, o livro proporciona diferentes olhares para a discussão na temática 
ambiental, com o objetivo proposto de reunir textos que procuram discutir o 
equilíbrio do meio físico. A obra ainda procura cumprir papel fundamental na 
disseminação de pesquisas nas geociências, principalmente para estudantes de 
Graduação e Pós-graduação. 

Dessa forma, convidamos a comunidade para conferir essa nova obra, 
com possibilidades e propostas para a aplicação e disseminação das pesquisas 
em diferentes contextos e realidades. 


Luis Ricardo Fernandes da Costa 
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RESUMO: A bacia hidrográfica do Rio das 
Velhas está localizada na região Central de 
Minas Gerais, com nascente no município de 
Ouro Preto, na Área de Proteção Ambiental 
da Cachoeira das Andorinhas. O Rio das 
Velhas é o maior afluente em extensão da 
bacia do Rio São Francisco, com 761 Km, 
da sua nascente até o ponto onde deságua 
no Rio São Francisco, em Barra do Guaicuí, 
distrito de Várzea da Palma (MG). Assim, o 
objetivo deste artigo é caracterizar a bacia 
hidrográfica do Rio das Velhas por seus 
principais aspectos físicos, populacionais e 
econômicos ao longo de todo seu território. 
Analisando na totalidade as características 
da bacia hidrográfica do Rio das Velhas, 
é notável a disparidade da distribuição 


Descobertas e desafios nas geociências 


Data de aceite: 02/05/2024 


da renda ao longo da bacia, retrato da 
diferença das atividades realizadas em 
seus diferentes territórios, tal qual a sua 
distribuição populacional. 

PALAVRAS-CHAVE: Análise 
Bacia hidrográfica. Minas Gerais 


integrada. 


ABSTRACT: The Rio das Velhas 
hydrographic basin is located in the Central 
region of Minas Gerais, with its source in the 
municipality of Ouro Preto, in the Cachoeira 
das Andorinhas Environmental Protection 
Area. The Rio das Velhas is the largest 
tributary in length of the São Francisco River 
basin, with 761 km, from its source to the 
point where it flows into the São Francisco 
River, in Barra do Guaicuí, district of Várzea 
da Palma (MG). Thus, the objective of this 
article is to characterize the Rio das Velhas 
hydrographic basin by its main physical, 
population and economic aspects throughout 
its territory. Analyzing the characteristics of 
the Rio das Velhas hydrographic basin in 
its entirety, the disparity in the distribution 
of income throughout the basin is notable, 
reflecting the difference in activities carried 
out in its different territories, as well as its 
population distribution. 
KEYWORDS: Integrated 
Hydrographic basin. Minas Gerais 


analysis. 
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INTRODUÇÃO 

A bacia hidrográfica do Rio das Velhas está localizada na região Central de Minas 
Gerais, com nascente no município de Ouro Preto, na Área de Proteção Ambiental da 
Cachoeira das Andorinhas. O Rio das Velhas é o maior afluente em extensão da bacia 
do Rio São Francisco, com 761 Km, da sua nascente até o ponto onde deságua no Rio 
São Francisco, em Barra do Guaicuí, distrito de Várzea da Palma (MG). Com área total de 
29.173 Km? (FEAM, 2002) abrange 51 municípios que concentram uma população de 4.8 
milhões de habitantes (IBGE, 2019) detentores de 62% do PIB do estado de Minas Gerais. 
O Rio das Velhas é a principal fonte de abastecimento de água na Região Metropolitana de 
Belo Horizonte e dos demais municípios que compõem a bacia, sendo essa subdividida em 
alto, médio e baixo curso (Guimarães, 1953). 

O nome “Rio das Velhas” foi dado pelo governador Antônio de Albuquerque Coelho 
de Carvalho no século XVII. Anteriormente o rio era chamado pelas tribos indígenas que 
ali habitavam de “Guaichui”, que significa “velhas tribos descendentes”. O fóssil humano 
mais antigo já encontrado no continente americano, conhecido como “Luzia” foi encontrado 
as margens do Rio das Velhas na região do município de Lagoa Santa. Segundo Goulart 
(2015), a bacia do Rio das Velhas se destaca no cenário mundial por abrigar resquícios 
arqueológicos do período da megafauna, datados de cerca de 12 mil anos atrás. 

Para Santos (1997), a consolidação de diferentes territórios é um processo histórico 
e dinâmico, que resulta na transformação da paisagem e no impacto ambiental. A paisagem, 
como síntese da produção do espaço geográfico, reflete as ações humanas passadas 
e presentes. As paisagens atuais são produtos de sua história e de seus processos de 
produção do espaço, e não são fixas, mas continuamente mutáveis. 

A bacia hidrográfica do alto Rio das Velhas, em Minas Gerais, é um exemplo de 
como o uso e ocupação do solo podem gerar conflitos ambientais. A região é um importante 
manancial de abastecimento para a Região Metropolitana de Belo Horizonte, mas também 
é alvo de atividades econômicas como mineração, agricultura e pecuária. Essas atividades, 
muitas vezes, são realizadas de forma insustentável, causando poluição e degradação 
ambiental. Os conflitos na bacia do alto Rio das Velhas são resultado da disputa entre 
diferentes interesses econômicos e sociais. De um lado, estão os setores que dependem 
da água para a sua atividade, como o abastecimento público e a geração de energia. De 
outro lado, estão os setores que exploram os recursos naturais da região, muitas vezes 
sem considerar os impactos ambientais (Lemos, 2019). 

Visto o posto, o objetivo deste artigo é: caracterizar a bacia hidrográfica do Rio das 
Velhas por seus principais aspectos físicos, populacionais e econômicos ao longo de todo 
seu território, compreendendo como se dá o uso de seus recursos hídricos através de sua 
extensão, desde a região de sua nascente no Quadrilátero Ferrífero, passando pela região 
metropolitana de Belo Horizonte, na região central do estado até atingir sua foz no Rio São 
Francisco no Norte de Minas. 
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METODOLOGIA 


Para a realização das análises foram levantados dados em alguns dos principais 
repositórios de dados disponíveis em termos de gestão da bacia e da gestão de meio 
ambiente do estado de Minas Gerais. O Comitê da Bacia Hidrográfica do Rio das Velhas 
mantém um repositório de dados de geoprocessamento de toda a bacia proveniente de 
anos de estudos desde sua criação, cnamado SIGA VELHAS. Deste repositório foram 
retirados os limites da bacia e as delimitações territoriais utilizadas para gestão da bacia 
hidrográfica. Outro repositório utilizado foi o IDE SISEMA da Secretaria Estadual de Meio 
Ambiente de Minas Gerais, onde foram levantados os principais dados referentes a área 
física da bacia, como clima, geologia, pedologia e vegetação. Os dados populacionais foram 
retirados do último censo demográfico do IBGE de 2022, porém, os dados socioeconômicos 
mais recentes utilizados são os do censo demográfico de 2010, pois dados atualizados 
ainda não foram disponibilizados nos portais de banco de dados do IBGE. Os dados foram 
processados e convertidos em produtos cartográficos através do software QGIS versão 
3.22, a fim de facilitar a leitura dos dados relevantes para a composição do artigo. 


ÁREA DE ESTUDO 


A bacia se localiza na região central do estado de Minas Gerais e se alonga no 
sentido norte a sul. Divide-se a bacia hidrográfica em trechos, sendo esses denominados 
“alto”, “médio” e “baixo”. O Alto Rio das Velhas se inicia no Quadrilátero Ferrífero, 
no município de Ouro Preto e estende seu limite até os municípios de Belo Horizonte, 
Contagem e Sabará. Parte do município de Caeté também faz parte da região, tendo seu 
limite definido pela serra da Piedade. O limite entre o Alto e o Médio Rio das Velhas é 
traçado pelo curso do Rio Paraúnas, maior afluente do Rio das Velhas, e segue até seu 
limite a sul, o Córrego Salobinho, limite entre os municípios de Curvelo e Corinto. Por fim, o 
Baixo Rio das Velhas compreende o município de Corinto, Gouveia, Presidente Kubitschek, 
Monjolos, Buenópolis, Joaquim Felício, Lassance, Várzea da Palma e Pirapora. É no limite 
entre esses dois municípios, em Barra do Guaicuí, no distrito de Várzea da Palma, onde o 
Rio das Velhas deságua no Rio São Francisco. 
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Mapa 1. Divisão da bacia hidrográfica do Rio das Velhas em trechos. 
Fonte: IGAM, 1998 


RESULTADOS E DISCUSSÕES 


Clima 


A classificação climática de Nimer (1972) leva em conta parâmetros quantitativos 
meteorológicos de precipitação de pluviosidade e temperatura para destacar os meses 
mais secos e mais úmidos de uma determinada região e classifica-la. Para tal, primeiro 
é considerado a temperatura média do mês mais frio, e posteriormente é considerado de 
acordo o grau de umidade, destacando a existência das estações mais secas e úmidas. 

A precipitação na bacia hidrográfica do Rio das Velhas decresce no percurso 
Norte a Sul, sendo os maiores valores presentes na região serrana de Ouro Preto, com 
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pluviosidade anual média registrada em torno dos 1800 milímetros. Ao longo da bacia, 
em sua porção oeste, nota-se uma depressão chuvosa e na porção oeste, ao longo da 
serra do Espinhaço, por influência topográfica, é registrado uma média anual na casa dos 
1700 milímetros anuais. Na região do Baixo Rio das Velhas, mais especificamente entre as 
cidades de Pirapora, Buenópolis e Várzea da Palma é registrado uma média pluviométrica 
em torno dos 1100 milímetros anuais. 

A temperatura média varia em torno de 18º no alto da bacia até 23º ao norte. Como 
exceção é notado a queda na média da temperatura nas proximidades dos municípios de 
Diamantina, ficando em torno de 19º e 20º. 

Portanto, o clima na bacia hidrográfica do Rio das Velhas de acordo com a 
classificação de Nimer (1972) se desenha da seguinte forma: 


600000 750000 
| A 
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—— Rio São Francisco 
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EB Úmido; de 3 meses secos, Subquente 
|| Úmido; de 3 meses secos, Mesotérmico Brando 


Semi Úmido; de 4 a 5 meses secos 


Semi ÚUmido; de 4 a 5 meses secos, Quente 


EEB semi Úmido; de 4 a 5 meses secos, Subquente 


Semi Úmido; de 6 meses secos, subquente 
EB Minas Gerais 
Brasil 


Classificação climática na bacia 
hidrográfica do rio das Velhas 
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Fonte: CBH Velhas, 2021. 
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Mapa 2. Classificação climática na bacia do Rio das Velhas segundo Nimer. 
Fonte: IGAM, 2021. 
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GEOMORFOLOGIA 


Segundo Moreira (2006), predomina na bacia do Rio das Velhas rochas das eras 
Arqueana e Proterozoica, divididas em três grandes grupos litológicos, sendo esses 
o Quadrilátero Ferrífero, o grupo Bambuí e a Serra do Espinhaço/Serra do Cabral. O 
Quadrilátero Ferrífero compreende o Alto Rio das Velhas na região dos municípios de Ouro 
Preto, Belo Horizonte e Santa Luzia, sendo delimitada pela Serra do Curral. O Quadrilátero 
Ferrífero é rico em rochas como o manganês, itabirito, quartzitos, filitos, ouro, e o minério 
de ferro. No grupo Bambuí, presente em grande parte do curso médio da bacia, estão 
presentes rochas sedimentares, carbonáticas e metapelíticas. Na porção baixa da bacia 
encontra se o grupo Espinhaço/Cabral, composto de quartzitos, quartzitos micácios, filitos 
e dolomitos. Destacam se nessa região as Serras do Espinhaço e do Cabral, por sua 
altimetria elevada em comparação as áreas mais rebaixadas da bacia (Moreira, 2006). 

De acordo com CPRM (2004) A bacia está inserida no bloco Brasília, sendo composta 
em sua maior extensão por gnáissicos granitoides em estado médio de metamorfismo em 
sua área exposta, visíveis nas regiões centro sul, meio leste e norte nordeste do estado 
de Minas Gerais, onde é possível visualizar litologias do Terciário e do Quaternário 
representando o período Cenozoico, sendo dividido segundo o CETEC (1983), em quatro 
unidades geomorfoló- gicas: Serra do Espinhaço, Planaltos do São Francisco, Quadrilátero 
Ferrífero e Depressão Sanfranciscana. 
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Mapa 3. Unidades geomorfológicas da bacia do Rio das Velhas. 
Fonte: CPRM, 2004. 


PEDOLOGIA 


Os solos da bacia do Rio das Velhas têm formação dada pela combinação de fatores 
que vão desde clima, embasamento litológico, biota e tempo. O embasamento litológico possui 
ligação direta com a morfologia do relevo da bacia, caracterizado pelas áreas mais aplainadas 
em sua porção central e sul onde se encontra o embasamento rochoso de predominância 
sedimentar, mais suscetíveis a processos intempéricos, enquanto na porção leste e sul estão 
as rochas mais resistentes ao intemperismo, abrigando assim as maiores elevações da região 
da bacia (Silva, 2019). De acordo com a classificação de solos da Embrapa (2005), as classes 
predominantes na bacia são Argissolo Vermelho Amarelo, Argissolo Vermelho, Cambissolo 
Háplico, Cambissolo, Latossolo Vermelho, Neossolo Flúvico, Neossolo Litólico e a presença de 
afloramentos rochosos na região do espinhaço, distribuídos da seguinte forma: 
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Mapa 5. Representação pedológica da bacia do Rio das Velhas. 
Fonte: EMBRAPA, 2005. 


VEGETAÇÃO 


De acordo com o CBH VELHAS (2004), restam poucas áreas de vegetação nativa 
ao longo da bacia do Rio das Velhas, principalmente no curso baixo e médio, onde a 
maior parte da vegetação nativa do cerrado, predominante ao longo da bacia, deu lugar a 
atividade agropecuária e agroindústria. As matas ciliares foram reduzidas a áreas restritas, 
o que agrava O assoreamento em pequenas bacias tributárias de importantes afluentes 
do Rio das Velhas. No geral, a crescente expansão das atividades agrárias suprime a 
vegetação ao longo da bacia, o que impacta diretamente no assoreamento, a quantidade 
de material orgânico nos corpos hídricos e na qualidade da água. Segundo o Ministério do 
Meio Ambiente (2006), a distribuição da vegetação ao longo da bacia hidrográfica pode ser 
representada através da seguinte composição: 
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Mapa 6. Representação da vegetação na área da bacia do Rio das Velhas. 
Fonte: MMA, 2005. 


DINÂMICA SOCIOECONÔMICA 


A dinâmica socioeconômica na bacia não se dá de forma homogênea, pela forma 
distinta na configuração dos municípios da região norte na foz do Rio, para região sul, 
principalmente na região metropolitana de Belo Horizonte. De acordo com Silva (2019), 
predomina no alto da bacia, na região metropolitana de Belo Horizonte as altas taxas de 
urbanização, devido a menor extensão territorial e o porte demográfico das cidades. Nas 
regiões médio e baixo os municípios possuem menor densidade demográfica e uma maior 


extensão territorial. 
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Segundo o CBH VELHAS (2021), a Região do Velhas possui 4.405.860 habitantes 
estando 4.291.475 em áreas urbanas e apenas 112.385 nas áreas rurais. Aproximadamente 
2.978.593 destes residem em Belo Horizonte e região metropolitana. Ao longo da região da 
bacia os municípios que a compõe são majoritariamente urbanos, sendo que a proporção 
da população urbana supera a rural na maioria dos casos, porém, nas regiões médio e baixo 
ainda se percebe uma significativa parcela da população que vive nas áreas rurais, tendo 
destaque para os municípios de Taquaraçú de Minas, Santana do Pirapama, Jequitibá e 
Presidente Juscelino. (Silva, 2019). 
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Mapa 7. Distribuição da população na área da bacia do Rio das Velhas. 
Fonte: IBGE, 2022. 
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Quanto ao desenvolvimento humano, de acordo com o Censo Demográfico de 2010 
do IBGE, todos os municípios que compõem a bacia hidrográfica estão minimamente no 
patamar médio, e alguns estão nos patamares alto e muito alto. Os maiores valores se 
encontram nos municípios de Belo Horizonte, com o IDHM muito alto (0,810) e Nova Lima 
(0,813). 
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Mapa 8. Indice de Desenvolvimento Humano por Municípios (IDHM) na área da bacia do Rio das 
Velhas. 


Fonte: IBGE, 2010.c 


A circulação financeira e o arranjo produtivo local se dão por meio do setor de 
serviços de enfoque tecnológico como a produção de softwares e biotecnologia na região 
da capital. Fora de Belo Horizonte destacam se a extração de pedra sabão na região de 
Ouro Preto e de joias e pedras preciosas na região de Nova Lima. Além disso, nessas 
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regiões destacam, de forma geral, a importância do setor de mineração como contribuição 
pra economia local, enquanto nas regiões do médio e baixo Rio das Velhas se destacam 
a produção agrícola, a pecuária, a indústria siderúrgica e o setor de serviços (AGÊNCIA 
MINAS, 2017). 


CONSIDERAÇÕES FINAIS 


Analisando na totalidade as características da bacia hidrográfica do Rio das Velhas, 
é notável a disparidade da distribuição da renda ao longo da bacia, retrato da diferença das 
atividades realizadas em seus diferentes territórios, tal qual a sua distribuição populacional. 
Na região do Alto Velhas, área onde se concentram as atividades de mineração e onde estão 
os grandes adensamentos urbanos na região da bacia estão também os maiores índices de 
desenvolvimento humano. Ao longo do médio alto, médio baixo e baixo Rio das Velhas se 
diminui o adensamento populacional e o índice de desenvolvimento urbano e as atividades 
econômicas se diversificam para atividades de agricultura e pecuária primordialmente. 
Nessa região estão os dois municípios de menor IDHM da bacia, Congonhas do Norte e 
Presidente Kubitschek, municípios na área da serra do Cipó, na cordilheira do Espinhaço, 
municípios de pequeno porte que tem o turismo como principal atividade econômica. 

É possível relacionar as características físicas à exploração dos recursos naturais 
da bacia desde o princípio da história de ocupação da região. Devido à exploração dos 
recursos minerais na região do Alto Velhas ali se concentraram os principais adensamentos 
populacionais ao longo da bacia desde o período colonial, primeiramente com Ouro Preto 
e anos depois com o surgimento da capital do estado e o desenvolvimento econômico e 
populacional que o surgimento de Belo Horizonte impulsionou nos seus arredores. Ao longo 
da bacia, nas regiões médias e baixa, devido ao aplainamento do terreno as principais 
atividades ali desenvolvidas estão relacionadas a agroindústria e a pecuária, com menores 
adensamentos urbanos que na porção alta da bacia. 
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RESUMO: Os tratamentos das biomassas 
pré e pós pirolise tem sido usados para 
melhorar as capacidades de adsorção do 
biocarvão. Os objetivos deste trabalho 
foram avaliar as alterações que a ativação 
do biocarvão de eucalipto por HNO, 
produziram nos componentes minerais, nos 
atributos físicos e químicos e na capacidade 
de adsorção dos elementos potencialmente 
tóxicos (EPT) Cd, Cr, Cu e Pb em soluções 
aquosas. Os experimentos testando a 
seletividade de adsorção dos mesmos EPT 
são raros. A ativação do BC com HNO, 
elevou a área superficial específica e a 
capacidade de troca de cátions e reduziu as 
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cinzas e a alcalinidade do BC, os: Mesmo 
com maior disponibilidade de grupos 
funcionais o BC, o; teve um desempenho 
inferior ao BC na adsorção dos ETP em 
experimentos com pH não tamponado. 
PALAVRAS-CHAVE: biocarvão de 
eucalipto, ativação com HNO,, adsorção, 
componentes minerais, elementos 
potencialmente tóxicos 


ACTIVATION OF EUCALYPTUS 
BIOCHAR FOR THE REMEDIATION 
OF SOLUTIONS CONTAINING Cd, Cr, 
Cu AND Pb 
ABSTRACT: Pre- and post-pyrolysis 
biomass treatments have been used to 
improve the adsorption capabilities of 
biochar. The objectives of this work were to 
evaluate the changes that the activation of 
eucalyptus biochar by HNO, produced in the 
mineral components, physical and chemical 
attributes and in the adsorption capacity of 
the potentially toxic elements (EPT) Cd, Cr, 
Cu and Pbin aqueous solutions. Experiments 
testing the adsorption selectivity of the same 
EPTs are rare. Activation of BC with HNO, 
increased the specific surface area and 
cation exchange capacity and reduced ash 


and alkalinity of BC, «os: Even with greater 
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availability of functional groups, BCHNOS performed worse than BC in the adsorption of ETP 
in experiments with unbuffered pH. 

KEYWORDS: eucalyptus biochar, activation with HNO3, adsorption, mineral components, 
potentially toxic elements 


INTRODUÇÃO 


O biocarvão possui um grande potencial para a remediação de águas e solos pelo 
seu baixo custo de produção, eficiência e resistência à degradação. As influências que as 
propriedades físicas, químicas, tipo de biomassa e tipo de pirolise tem sobre as capacidades 
de adsorção de metais e ânions pelos biocarvões têm sido intensamente pesquisadas (Tan 
et al., 2015, Liu et al. 2022; Nguyen et al., 2023). Além dos modos de produção vários 
tratamentos na biomassa pré e pós pirólise tem sido realizados com o intuito de melhorar 
sua capacidade de remediação de solos e águas. Os métodos mais comumente utilizados 
para modificar e aumentar a eficiência do biocarvão na remoção de poluentes do solo e 
da água são: tratamento ácido ou oxidação (remoção de matéria orgânica não pirolizada 
e das cinzas dos poros, aumento da área superficial e dos grupos funcionais carboxyl) e 
básico (remoção e matéria orgânica não pirolizada e das cinzas dos poros, aumento da 
área superficial e dos grupos funcionais hydroxyl), amination, ativação por vapor (remoção 
e matéria orgânica não pirolizada dos poros com aumento da área superficial e dos grupos 
funcionais hydroxyl), produção de biocarvão com cobertura de óxidos (aumento das cargas 
superficiais positivas), matéria carbonácea, argilas (aumento das capacidades de adsorção 
de oxyanions e cátions poliatômicos), compostos orgânicos (aumento da capacidade e da 
energia na adsorção de metais) e biofilmes de grafeno e nanotubos de carbono (aumento 
da área superficial e dos grupos funcionais contendo oxigênio) (Rajapaksha et al., 2016; 
Sizmur et al., 2017; Lietal., 2017; Huang et al., 2021, Hamid et al., 2022; Wang et al., 2022; 
Yu et al., 2022). 

Biocarvões ativados com ácidos tem sido utilizados com sucesso na adsorção de 
fons em soluções. O biocarvão de fibras de cactos O. ficus indica, ativado com HNO, 
removeu Cu em soluções com pH - 6, por meio da adsorção pela formação de complexos 
de esfera externa (Hadijittofi et al.,2014). Esse mesmo biocarvão removeu com grande 
eficiência o U(VI) tanto em soluções ácidas (pH 3) como neutras (pH 6.5) por meio tanto da 
formação de complexos de esfera interna como de esfera externa (Hadiittofi & Pashalidis, 
2015). O biocarvão de palha de milho, ativado por diferentes ácidos (HNO,, H,SO, e 
H.;PO,), teve uma elevação na capacidade de remoção de Cr (VI) quando comparado com 
o biocarvão in natura. A ativação ácida aumentou os grupos funcionais fenol e hydroxyl, o 
que favoreceu a redução do Cr(VI) em Cr (Ill) uma vez que estes grupos funcionais que 
se comportaram como doadores de elétrons. O biocarvão ativado por HNO, apresentou a 
melhor capacidade de remoção entre os pH's 3-9. O processo de redução do Cr foi mais 
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eficiente em baixos pH's (Xu et al. 2020). Zhang et al. (2020) observaram um aumento no 
sequestro de Pb em biocarvões de lodo de esgoto tratados com HCl e HF em comparação 
com o biocarvão in natura. 

Apesar dos casos de sucesso citados, a remoção das cinzas pelos tratamentos de 
ativação ácida pode provocar efeitos antagônicos na capacidade de remoção de íons em 
solução pelos biocarvões. A fração mineral do biocarvão, presente em suas cinzas, como 
os metais alcalinos e alcalinos terrosos em carbonatos, fosfatos ou óxidos solúveis em 
água, tem influência na capacidade de adsorção dos biocarvões e seu uso (Xu et al., 2017; 
Wang et al., 2015; Meng et al., 2022; Ye et al., 2022), podendo ser responsável por até 
99% da adsorção de Pb (Xu et al., 2017) e favorecer a remoção de Zn (Qian et al., 2016). A 
remoção de Cu de soluções por biocarvão de casca de amendoim com elevados teores de 
cinzas, que continham PO4 3-, COS 2-, OH- e SO4 2-, foi mais eficiente do que o biocarvão 
que teve as cinzas removidas, tendo havido provavelmente a precipitação de sais de Cu 
com estes ânions, em especial o PO43- (Zhou et al., 2016). A adsorção e a precipitação de 
possíveis sais com metais pesados, a partir da solubilização das cinzas dos biocarvões, 
pode ser investigada por meio da modelagem geoquímica, onde é possível identificar os 
principais sítios ou grupos funcionais no processo de adosrção (Kypritidou et al., 2022) e a 
saturação e precipitação de sais (Parkhurst & Appelo, 2013). 

Uma vez que a ativação química ácida e a remoção das cinzas podem produzir 
resultados antagônicos nas propriedades e capacidades de retenção de metais dos 
biocarvões, neste trabalho foram avaliadas as alterações que a ativação do biocarvão 
de eucalipto por HNOS produziram: (i) na alcalinidade e nos componentes minerais do 
biocarvão, (ii) nos atributos físicos e químicos e (iii) na capacidade de adsorção dos 
elementos potencialmente tóxicos Cd, Cr, Cu e Pb em soluções aquosas. 


MATERIAIS E MÉTODOS 


O biocarvão de eucalipto (BC) foi produzido com cascas de eucalipto, com 
tamanho entre 0.1-0.5 mm, por pirolise rápida (3s) em temperaturas entre 450-500ºC, 
foi obtido 
a partir do tratamento químico do BC: 50g do BC foi colocado num béquer com 1L de 


em condição aeróbica. O biocarvão de eucalipto ativado com HNO, (BC, os) 
solução 1 molar de HNO, (HNO, 65% P.A.). A mistura foi mantida em agitação a 80 ºC 


por 2 horas em agitador magnético. Posteriormente o BC foi filtrado a vácuo para a 


HNO3 
retirada da solução remanescente e lavado com água desionizada (Millipore?) até não 
ser mais possível identificar o NO, nitrato na água de lavagem (determinação de N por 
destilação em Destilador Kjeldahl e titulação). Por fim foram secas em estufa a 60ºC até 
peso constante por + 2 dias. 

No BC e no BC,,,o; foram determinados o pH em água (1:10); a umidade a 65ºC; 


os teores de cinzas (800ºC — ASTM), o C fixo e os voláteis. Os teores totais de C, Ne H 
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foram determinados num analisador elementar CHN Perkin Elmer 2000, e o O foi obtido por 
diferença (0=100-C-H-N). Foi realizada solubilização parcial das cinzas segundo o método 
EPA3051a (USEPA, 1998) com leitura em ICP-OES para a determinação dos teores 
semi-totais de macro e micronutrinetes metálicos e elementos potencialmente tóxicos. 
A capacidade de troca de cátions (CTC) e o ponto de carga zero (PCZ) foram obtidos 
segundo descrito em Lee et al.(2010) e UCHIMIYA et al. 2011, respectivamente. O BC e 
no BC 
Phaser ; (ii) por microscopia eletrônica de varredura (MEV) no microscópio eletrônico LEO 


unos foram analisados por: (i) difração de raios X (DRX) no equipamento Bruker-D2 
430i-Oxford e (iii) por espectrometria no infravermelho com transformada de Fourrier (FTIR) 
num Agilent Cary 630 FTIR spectrometer. As áreas superficiais foram determinadas pelo 
método Brunauer-Emmett-Teller (BET — N,). 

Os experimentos de adsorção foram realizados adicionando-se em tubos de centrífuga 
de 50 mL 0,25 g de adsorvente (BC ou BC 
Cu ou Pb) em condições de pH 4 e pH 5, contendo os cátions nas concentrações de 0,1 


HNO) e 40 mL de soluções monometálicas (Cd, Cr, 
a 1,8 mM. O intervalo de pH dos experimentos foi escolhido com base na capacidade de 
solubilização dos metais e por ser encontrado em águas subterrâneas afetadas por chuvas 
ácidas (Kjoller et al., 2004) e rios poluídos (Oberholster et al., 2017). As soluções foram 
preparadas a partir de diluições soluções monometálicas 10 mM dos reagentes (PA) Cd 
(NO,).2H,0; CuCl,. 2H,0; CrCl,. GH,O0 e Pb (NO,),. 6H,0, usando-se água desionizada. Os 
experimentos foram executados com amostras em duplicata. As amostras foram agitadas 
em agitador na vertical, à 210 oscilações por minuto por 24 horas a 24ºC. Ao término da 
agitação as amostras foram centrifugadas e filtradas com filtro de acetato de celulose de 
0,45 um. Os extratos foram preservados com ácido nítrico e analisados num ICP-OES MPX 
Varian para a determinação dos íons em solução. Foram utilizados brancos, duplicatas 
e padrões internos do laboratório como controle de qualidade analítico. A descrição dos 
experimentos de adsorção pelos modelos de Langmuir e Freundlich utilizaram as equações 
apresentadas em Limousin et al. (2007). 

O software PHREEQC (Parkhurst e Appelo, 2013) foi utilizado para realizar 
modelagem geoquímica de soluções, para cálculos de força iônica, determinação das 
especiações de seus íons e cálculo de índices de saturação. 


RESULTADOS E DISCUSSÕES 


A pirólise rápida das cascas de eucalipto produziu um biocarvão (BC) com teores 
CHNO de 60%, 3,8%, 1% e 35% respectivamente, com 12,5% de cinzas e rendimento 
em torno de 20% em relação a biomassa inicial (TABELA. 1). A ativação do BC após sua 
produção com o HNO,, provocou a redução do H (3.4%), aumentou o teor de N (2,5%), 
manteve inalterado o teor de C (60%) e O (35%) do BC 


HNO3* 
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Biocarvão BC BC- HNO, 

C% m/m 60 +1 60 +1 

H% m/m 3,87 + 0,09 3,4 +0,1 

N% m/m 1,04 + 0,04 2,54 +0,05 

O % m/m 35,09 34,06 

H/C 0,774 0,68 

O/C 0,44 0,43 

N/C 0,015 0,04 

Fórmula Empírica CH, 570,47 Noo7 Css so O, Nous 

ASE m?g 2,759 7,583 

CTC (mmolc/kg) 445 604 

pH 7,9 3,2 

pH equil. c/ Sol pH 4 7,2 3,6 

pH equil. c/ Sol pH 5 7,2 3,8 

PZC 2,04 2,03 

Umidade % m/m 6,6+0,1 7,6+0,1 

Cinzas % m/m 12,4+0,5 6,3+0,1 

Material volátil % m/m 42,7+0,1 52,4+0,1 

Carbono fixo % m/m 45,8+0,8 41,3+0,2 
Extração Parcial — Método EPA 3051 

mg/kg BC BC, nos 

Al 367,5+6,27 406,5+31,1 

As 4,3+0,75 0,1+0,1 

Ba 64,6+2,61 37,6+0,6 

B 16,9+0,40 0,00,0 

Ca 6400+0,20 1000,0+0,1 

Pb 7,4+0,05 2,5+0,4 

Cu 10,3+0,05 3,7+0,3 

Cr 27,4+0,09 10,2+1,2 

s 800+0,00 200+0,0 

Fe 20,8+1,30 977,8+22,8 

P 2700+0,03 300+0,1 

Mg 3000+0,04 400+0,0 

Mn 121,0+0,13 64,7+4,1 

Ni 8,5+0,35 1,1+0,5 

Se 7,0+3,25 2,5+0,9 

Zn 58,6+3,61 21,5+3,5 

K 41187,8+614,06 328,2+0,2 

Na 288,3+0,06 15,8+0,2 


Tabela 1. Atributos do físicos, físico-químicos e químicos BC e BC 
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As áreas superficiais do BC e BC apresentaram valores de 2.7 m?g e 7.6 m?/g, 


HNO3 
respectivamente. As diferenças das ASE se manifestaram na CTC: o BC apresentou CTC 


de 445 mmolc kg! enquanto no BC a CTC aumentou para 604 mmolc kg'!. Como o 


HNO3 
tratamento do BC com o HNO, diminuiu os teores de cinzas de 12% para 6% houve também 


a diminuição do pH, de cerca de pH 8 no BC para pH 3 no BC 


HNO3” 
O tratamento do BC com HNO, teve como objetivo de promover a oxidação do 


BC e aumentar a polaridade e a disponibilidade dos grupos funcionais na superfície do 
BC, por meio da saturação das ligações do C e inserção de heteroatomos (H,O,N) nestes 


grupos funcionais. As razões O/C (Tabela 1) indicam que no BC não houve o aumento 


HNOo3 
esperado da polaridade decorrente da ligação de O aos C's. Apesar do pequeno aumento 


da polaridade no BC houve o aumento da CTC provavelmente devido a elevação da 


HNO3º 


ASE. O teor dos voláteis aumentou no BC provavelmente em função destruição cadeias 


HNO3 
de C pelo ataque ácido, coerente com a relação molar H/C que diminuiu de 0,76 para 0,68, 


o que geralmente está associado a um maior caráter aromático do BC os: 


EXTRAÇÃO PARCIAL DAS CINZAS- MÉTODO EPA 3051A 

Na extração parcial das cinzas do BC obtida pelo método EPA3051a, os íons com 
maiores concentrações foram K, Ca, Mg, P e em menores teores S, Al, As, Pb, Cu, Cre 
Ni (TABELA 1). No BC 
e a redução dos teores de cinzas (TABELA 1). Em virtude deste menor teor de cinzas, 


pnos à sua ativação acida com HNO, promoveu a dissolução parcial 
houve a diminuição dos teores de K ,Ca, Mg, Na, Pe S obtidos na extração parcial com 
o método EPA 3051a, uma vez que estes elementos ocorrem no conteúdo mineral das 
cinzas. Na ativação houve a troca dos cátions presentes nos complexos de adsorção do 
BC 
também a diminuição das concentrações de outros íons extraídos como As, Pb, Cu, Cr, e 


unos Pelo H*, o que deve ter contribuído para diminuição do seu pH (TABELA. 1). Houve 
Ni, possivelmente associados tanto às cinzas quanto aos sítios de troca. 

Segundo Fidel et al. (2017) a alcalinidade total dos biocarvões se correlaciona 
fortemente com as concentrações de cátions básicos. Entretanto a alcalinidade total é uma 
função dos grupos funcionais (estrutural), dos carbonatos, dos álcalis orgânicos solúveis e 
de outros álcalis inorgânicos, e vária também segundo a matéria prima e modo de produção. 
A ativação do BC com HNO, reduziu tanto o teor de cinzas como o teor dos cátions básicos 
(TABELA 1), o que resultou na redução de sua alcalinidade total, que se manifestou no pH 
do BC, oq: de pH 7,9 (BC) para pH 3,2 (BC 


HNO3" HNo3)- 
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CARACTERIZAÇÃO DOS BIOCARVÕES POR DRX, MEV E FTIR 


No BC foram identificados por DRX sólidos cristalinos com a estrutura semelhante 
à do quartzo e da silvita no BC (Fig 1). No BC 
estrutura do SiO, (Fig 1). O KCl identificado no BC deve fazer parte das cinzas pois não foi 


unos foi encontrado apenas o cristal com 


encontrado no BC, o: que teve seu teor de cinzas diminuído, provavelmente tendo sido 
eliminado pelo tratamento com os ácidos. No BC são observadas estruturas celulares de 


tecidos vegetais preservadas (Fig 2) enquanto que nos BC as partículas são menores e 


HNO3 
as estruturas menos visíveis, em virtude do ataque ácido sofrido no tratamento. 


— BE UC 
— BEUC-HNO3 


Angulo 20 


Figura. 1 Difração de raios X do BC e do BC, os- 

No exame por FTIR o interferograma do BC (Fig. 3) apresentou uma assinatura 
compatível com a da lignina (Guo et al., 2008), presente nas cascas de eucalipto in natura. 
No BC (Fig S3 e TABELA. S1) foram identificados picos relacionados com estiramentos 
de ligações O-H (grupos hidroxila associado a compostos fenólicos relacionados com a 
lignina), C-H (alifático), C=O (grupo carboxíla), C=C (aromático, indicativo de lignina), 
C-O-C (celulose), C-O (fenólico). No BC, os 
em 1426 e o aparecimento das bandas de absorção C=O em 1710 e C-O em 1261. A 


há a ausência da banda de absorção C=C 


eliminação dos picos de C=C no BC sugerem a oxidação do C e a diminuição dos 


HNOSº 
anéis aromáticos, ou pelo menos das ligações saturadas do C. Além disto, a intensidade do 
pico C=C em -1600 diminuiu no BC, o, possivelmente devido à saturação de parte destas 


ligações. 
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bg 1 2 3 [4 |5 
o | 58 74 
si | 40 11 
ci |0,42 8.2 
Kk [114 16.6 
Mg - - 

Total | 100 100 


B Cc 


Figura. 2: A— BC Imagens de MEV mostrando estruturas celulares preservadas de tecidos de 
plantas, ornamentos de SiO2 e pontos de análise. B and C — Imagens de MEV do BC, o; Mostrando 
ornamentos de SiO, e pequenas particulas e estruturas parcialmente destruidas devido ao ataque acido 
sofrido no tratamento de ativação. 
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Guo et al. (2008) 
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Figura 3 — Interferogramas do BC e do BC, o; (este trabalho) e da lignina de Guo et al. (2008). 


ISOTERMAS DE ADSORÇÃO 


Os resultados dos experimentais foram bem descritos pelos modelos de Langmuir e 
Freundlich (r,? e r?>0,80) na maioria dos experimentos (TABELA 2 e Fig. 4). As capacidades 
máximas de adsorção obtidas pelo modelo de Langmuir do Cd, Cu e Pb no BC em pH4 
variaram entre 0,30-0,40 mmol g”, enquanto a do Cr foi de 0,29 mmol g (TABELA 2 e Fig. 1). 
A ordem decrescente de capacidade máxima de adsorção encontrada foi: Cd>Cu>Pb>Cr. 
Na solução de metais com pH 5 houve o aumento de qmax do Cr (0,77 mmol g) e do 
Pb (0,47 mmol g!) enquanto o Cd e o Cu apresentaram diminuições em seus valores de 
qmax, 0,23 e 0,36 mmol g'! respectivamente. Desta forma, em pH 5 a ordem de qmax é 
Cr>Pb>Cu>Cd. Em virtude da provável participação das cinzas na retenção dos íons via 
precipitação, não é possível atribuir significado para qmax e n. 
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Langmuir 


qmax (mmol/g) K, (L/9) re 
pH4 pH5 pH4 pH5 pH4 pH5 
Cd 0.42 0.18 15.8 91.6 0.96 0.97 
Sé Cr 0.29 0.94 4.3 4.8 0.98 0.99 
Cu 0.40 0.26 11.2 39.4 0.88 0.82 
Pb 0.30 0.45 579 3.6 0.47 0.82 
Cd 0.11 0.09 6.2 4.3 0.88 0.83 
Cr 0.01 0.52 105 3.7 0.86 0.97 
BC, nos 
Cu 0.13 0.16 7.5 4.2 0.94 0.91 
Pb 0.11 0.12 90.8 106 0.93 0.81 


r,?= coeficiente de correlação; qmax = capacidade máxima de adsorção; K, = entalpia relacionada 
com a constante de adsorção. 


Freundlich 
K(L/9) :s 
pH4 pH5 pH4 pH5 pH4 pH5 
Cd 0.41 0.23 0.58 0.24 0.95 0.88 
BC Cr 0.44 0.77 0.25 1.85 0.95 0.99 
Cu 0.44 0.36 0.55 0.44 0.96 0.93 
Pb 0.12 0.47 0.43 0.43 0.54 0.81 
Cd 0.38 0.35 0.10 0.07 0.93 0.93 
Cr 0.16 0.59 0.01 0.55 0.51 0.98 
BCnos 
Cu 0.37 0.40 0.11 0.12 0.92 0.94 
Pb 0.19 0.24 0.13 0.15 0.84 0.95 


rp?= coeficiente de correlação; n = diversidade de sítios de adsorção; K, = constante de capacidade 


de sorção 


Tabela 2. Parâmetros ajustados de experimentos de adsorção de Cd, Cr, Cu, Pb em soluções 
monometálicas com pH 4 e 5 pelos BC e BC, o; pelos modelos de Langmuir e Freundlich. 


Aalteração da intensidade da adsorção do Cr entre os pH's 4 e 5 foi a maior alteração 


observada na ordem de adsorção. A diversidade dos sítios de adsorção, representada 


pelos grupos funcionais identificados no FTIR (Fig. 3), também se manifestou por meio do 


parâmetro n modelo de Freundlich (TABELA 2) que apresentou valores inferiores a 0,50. 
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Figura 4. Isotermas de adsorção de Langmuir e Freundlich para os experimentos de adsorção de 
soluções monometálicas dos EPT"s pelos BC e BC, os: 
previamente aos experimentos de adsorção eram 4 e 5. 


As capacidades de adsorções máximas dos EPT pelo BC e pelo BC 


Os pH's das soluções monometálicas 


pnos estão dentro 


do intervalo encontrado para os biocarvões provenientes de diversos tipos de biomassas 
produzidos em diversas temperaturas (TABELA 3): Cd de 0,3 a 51 mg L?, Crde 3 a 123 mg 
L!, Cu de 4,6 a 54 mg Le Pb de 2 a 147 mg L'. A ordem de adsorção se altera quando se 
usa a massa no lugar da molaridade (TABELA 2 e TABELA 3): BC-pH4 — Pb>Cd>Cu>Cr, 
BC, o; PH4 — Pb>Cd>Cr>Cu, BC-pH5 — Pb>Cr>Cd>Cu e BC, o; PH5 — Cr>Pb>Cd=Cu. 
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Capacidade de adsorção (mg/9) 
T(ºC) Cr Cd Pb Cu 
Lignina — material sem queima - - 29,7 89,4 22,6 Guo et al., 2008 
BC — pH4 450-500 is! 47 62,1 25,2 Este trabalho 
BC, nos PH4 450-500 |0,52 |123 |227 |82 | Estetrabalho 
BC — pH5 450-500 | 49 20,2 98,2 16,3 Este trabalho 
BC,nos” PHS 450-500 | 24 10,1 24,9 10,1 Este trabalho 
Lascas de madeira 200 25,5 39 147 52 El-Shafey et al., 
2002 
Palha de gergelim 700 65 86 102 55 Park et. al., 2016 


Tabela 3: Capacidade máxima de adsorção dos EPT's no BC, no BC, o;e em outros biocarvões 


Nos raros experimentos testando a seletividade de adsorção dos mesmos EPT 
deste trabalho, a partir de soluções aquosas para um mesmo biocarvão, foram encontradas 
as seguintes sequencias de adsorção (TABELA 3): Pb>Cu>Cd>Cr (El-Shafey et al., 2002) 
e Pb>Cd>Cr>Cu (Park et. al., 2016). Para a lignina, polímero precursor das paredes das 
células das plantas terrestres e presente na casca de eucalipto (Andrade et al., 2010) com 
assinatura semelhante à e BC no interferograma do (Fig. 3), Guo et al. (2008) encontraram 
a seguinte sequencia de adsorção: Pb>Cd>Cu, que foi semelhante a todos os experimentos 
aqui realizados com exceção do BC, o, PH5. Os valores máximos de adsorção da lignina 


são muito próximos aos encontrados para o BC-pH5 (TABELA. 3) e superiores ao BC, os: 


MECANISMOS DE ADSORÇÃO 


A modelagem hidrogeoquímica indicou que nas soluções de adsorção, nos pH's 4 
e 5, as espécies de Cd, Cu e Pb ocorrem na forma de íons positivos bivalentes hidratados. 
Nos pH's de equilíbrio destas soluções com os biocarvões a especiação destes íons 
continuou a mesma (TABELA. 4). O Cr por sua vez ocorreu como o íon Cr*, que devido 
a sua elevada eletronegatividade provoca a hidrolise das moléculas de água formando 
oxiânions, nas condições deste trabalho o íon predominante foi o HCrO. 
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Cd Cr Cu Pb 
Solução Cd? HCrO Cu*>>>CuCl*>CuOH* | Pb'2>>PbNO,*>>PbOH* 
monometálica - >>Cr,0,*> 
pH4 CrO,? 
Solução Cdt HCrO >> | Cut>>> CuOH* = CuCl* | Pb?>>PbNO,*>PbOH* 
monometálica - Cr,O,?= 
pH5 Cro,? 
Solução em Cd? HCrO >> Cu*>>>CuCI* Pb'>>PbNO,* 
equilíbrio com Cr,O? 
BC, nos. PH-3,5 
Solução em Cd? HCrO; Cu>> CuOH*> CuCO, | Pb*>PbNO,'>PbOH* 
equilíbrio com BC =CrO,*> 
- pH-6.5 Cr,O? 


Sólidos com índices de saturação > 1 nas soluções em equilíbrio com os biocarvões 


Solução em 

equilíbrio com Todos insaturados 

BC, nos. PH-3,5 

Solução em Cd,(PO4), nenhum Cu4CI2(0H),, Cu,(PO4),, | Pb.(PO,).CI, Pb5(OH) 


equilíbrio com BC 
- pH-6.5 


Cu,(PO4),:3H20, 
Cu,CO,(OH),, CuO 


(PO4),, Pb.(PO ),, 
PbHPO,, PbCO,, 


PbO,, Pb3(COS) (OH), 
Pb4CI2(OH), 


Tabela 4. Espécies iônicas nas soluções de adsorção antes e após os experimentos de adsorção com 
oBCeoBC e indices de saturação (IS) de sólidos nas soluções de adsorção em equibilibrio com 
oBCeoBC 


HNO3 * 
HINOS * 


Em virtude da similaridade dos íons bivalentes em solução, os mecanismos de 
adsorção pelos biocarvões seriam compartilhados entre eles (Li et al., 2017; Sizmur et 
al., 2017): a troca iônica estequiométrica, a atração eletrostática, a complexação e a 
precipitação na superfície. A troca iônica e a complexação são dependentes do pH do meio 
e ocorrem nos grupos funcionais, como o fenólico, carboxílico e hidroxílico, encontrados 
no BC e BC,,o; (TABELA 2, Fig.3). A adsorção por atração eletrostática se dá entre íons 
com cargas positivas e nuvens eletrônicas deslocadas associadas a grupos aromáticos 
criando interações cátion-m com as ligações aromáticas C=C (Li et al., 2017). A precipitação 
dos íons na superfície do biocarvão ocorre em virtude da formação de sais com baixa 
solubilidade formados a partir dos anions liberados pelas cinzas do biocarvão, como PO *, 
CO,”, SO? e os EPT's. 

Dentre os EPT aqui estudados o Cr apresenta um mecanismo de adsorção 
diferenciado. Os modelos de adsorção do Cr** em espécies negativamente carregadas, 
como HCrO,, indicam (i) a atração eletrostática entre estas espécies e as superfícies 
carregadas positivamente dos biocarvões e (ii) a redução do Cr** para o Cr** principalmente 
nos grupos funcionais que contém oxigênio, como o carboxílico e o fenólico, e posterior 
complexação do Cr** pelos grupos funcionais (Li et al., 2017). 

A solubilização das cinzas do BC elevou os pH's de equilíbrio das soluções de 
adsorção em contato com o BC para valores em torno de pH 7,2. Nestas condições há 
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o favorecimento da desprotonação dos grupos funcionais, aumentando a quantidade de 
cargas negativas do BC e aumentando a capacidade de adsorção de cátions. Além disto, os 
íons liberados pela dissolução das cinzas (K, Ca, Mg, Cl, P) também podem ter favorecido 
a precipitação de sais insolúveis. Na modelagem hidrogeoquímica vários compostos de 
Cd, Cu e Pb com os possíveis íons liberados pela dissolução das cinzas apresentaram 
índices de saturação maiores do que 1 (Tabela 4) para a solução com pH=7, o que indica a 
possibilidade de ter havido a precipitação destes compostos. 


Por outro lado, no BC o tratamento recebido com o HNO, para a sua ativação 


HNO3 
promoveu a solubilização parcial das cinzas e a remoção dos seus íons dissolvidos. Como 


consequência o pH de equilíbrio do BC com as soluções contendo os metais dos 


HNo3 


experimentos de adsorção ficou entre 3,6-3,8. Como o PCZ do BC é pH 2,08, na sua 


HNO3 
superfície já há existência de cargas negativas, possibilitando a adsorção de íons positivos. 
Porém, os principais sítios presentes para a adsorção dos EPT na lignina e encontrados 
no BC e BC,os 


favorecidas pelo aumento do pH. Este comportamento foi observado por Guo et al. (2008) 


o fenólico e o carboxílico, possuem capacidades adsortivas que são 


em seus experimentos, onde a diferença na capacidade de adsorção entre os pH's 3,5- 
5,5 foi superior a 2 vezes para Cu e Pb e 8 vezes para o Cd. Em virtude do baixo pH de 
equilíbrio do BC 
formadas na superfície do BC 


unos COM as soluções dos experimentos, a quantidade de cargas negativas 


unos deve ter sido pequena e limitou sua capacidade de 


adsorção. Além disto, como não há a disponibilidade dos íons liberados pela dissolução 


das cinzas, a precipitação de sais insolúveis dos EPT no BC não é esperada. Estes dois 


HNO3 


fatores fizeram com que o BC apresentasse uma capacidade de adsorção menor que o 


HNO3 
BC, mesmo possuindo uma CTC maior. 


CONCLUSÕES 


A ativação do BC com HNO, elevou a ASE e a CTC e reduziu as cinzas e a 
alcalinidade do BC, os 
teve um desempenho inferior ao BC na adsorção dos ETP em experimentos com pH não 


Mesmo com maior disponibilidade de grupos funcionais o BC, os 
tamponado. O teor de cinzas e os outros fatores que influenciam na alcalinidade total 
do BC devem ter tido uma maior importância na adsorção dos EPT, em virtude do maior 
pH de equilíbrio das soluções, favorecendo a formação de cargas superficiais negativas 
em seus grupos funcionais e a precipitação de sais insolúveis dos ETP. Desta forma, a 
queda o desempenho do BC em sistemas de remediação de soluções aquosas contendo 
metais, em virtude principalmente do consumo de sua reserva alcalina, deve ser levada em 
consideração na concepção de sistemas de remediação. 
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LIMITED HEIGHTS OF VERTICAL CLIFFS AND 
MOUNTAIN WALLS LINKED TO FRACTURING IN DEEP 
TUNNELS - Q-SLOPE APPLICATION IF JOINTED SLÓOPES 
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ABSTRACT: Intact brittle rock can fail in 
tension even when all principal stresses 
are compressive. This is due to lateral 
expansion and extension strain when near 
to a free surface, caused by Poisson's ratio. 
Exceeding tensile strength due to stress 
anisotropy and Poisson's ratio are the 
fracture-initiating conditions around deep 
tunnels, not the increasing mobilization of 
compressive strength, commonly beyond 
0.4x UCS. In a related discovery, the limiting 
height of vertical cliffs and near-vertical 
mountain walls can also be explained using 
extension strain theory. The range of limiting 
heights of approximately 20m for cliffs in 
porous tuff to record 1,300m high mountain 
walls in granite are thereby explained. Tensile 
strength is the weakest link behind cliffs and 
ultra-steep mountain walls. Sheeting joints 
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can also be explained by extension strain 
theory. Maximum shear strength is the 
weakest link when stress levels are ultra- 
high, or when there is jointing and maximum 
slope angles is the issue. Here one can use 
Q-slope. The world's highest mountains are 
limited to 8 to 9km. This is due to non-linear 
critical state rock mechanics. It is not due 


to UCS. 
KEYWORDS: Deep  tunnels, Cliffs, 
Mountains; Extension strain; Tensile 


strength; Shear strength 


INTRODUCTION 


The lessons from fracturing in 
deep tunnels is the starting point for the 
ultra-simple cliff-height and mountain 
wall-height equation which is introduced 
in this article. The observed and recently 
modelled fracturing behavior of deep 
tunnels in massive rock indicates that 
fracturing may be initiated by extensional 
strain over-coming the tensile limit, even 
when all stresses are compressive. This 
is possible due to the lateral expansion 
caused by Poisson's ratio. A small-scale 


example of this is the acoustic emission 
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that occurs due to micro-fracture initiation when testing intact rock cylinders in traditional 
uniaxial compression, where Poisson's ratio is also at work. The commontly used parameter 
obtained from such tests is 0,, the unconfined compression strength (commonly written as 
UCS). This might be 150MPa for granite but only 1.5MPa for weak porous tuff, the medium 
once used by Christian cliff-dwellers in Cappadocia, Turkey. The tuffs are so weak that there 
have been many historic cliff failures, which expose old dwellings and Christian churches at 
irregular intervals. The most basic strength parameter o, has traditionally been compared 
with the estimated maximum tangential ('arching”) stress, to investigate if a deep tunnel will 
suffer fracturing or rock-burst and need more support like sprayed concrete and rock bolts. A 
newly excavated tunnel results in a big contrast between the maximum tangential ('arching”) 
stress (0,) and the almost unloaded radial stress (0,). For elastic isotropic materials and a 
circular tunnel, the theoretical maximum tangential stress is three times the major principal 
stress (0,) minus the minimum principal stress (0,) acting in the same plane, at right angles 
to the tunnel. At 1,000m depth we might have 0, = 3 x 30 - 15 = 75MPa, due to assumed 
in situ stresses o, (vertical) and o, (horizontal) of 30 and 15 MPa. If we now compare the 
magnitude of o, with the available uniaxial strength 0,, say 150MPa for granite, then the 
ratio 0,/0,= 75/150 = 0.5, will suggest from Figure 1 (and from Table 1) that fracturing and 
break-out may occur: i.e. more support is needed due to much increased SRF i.e. lower Q. 
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2, +aRC 
O Kirsten & Klokow, 1979 
O Jiayou et al, 1991 
1.8! O Martinetal, 1994 
W Martin, 1989 
O Ortepp & Gay, 1984 
16 X Pelli et al, 1991 
À Stacey & de Jongh, 1977 


0.49(+0.1) + 1.25% 


Figure 1. The traditional assumption is that fracturing in a deep tunnel is limited if compressive strength 

is sufficient. Martin et al. 1998, and Grimstad and Barton, 1993 have independently confirmed that when 
the stress/strength ratio 0,/0, > 0.4+0.1, fracturing and break-out is likely, as illustrated here, to a depth 

of 2 to 3m. This ocurred in several large diversion tunnels in Brazil (Ita HEP, Barton and Infanti, 2004). 


In the Q-system ratings for the stress reduction factor SRF, the rating is accelerated 
when passing a ratio 0,/0, of 0.4, as shown in Table 1. The numerous case records collected 


by Grimstad for the case of deep road tunnels in Norway, giving stress/strength ratios in 


the case of tunnels of up to 1.4km depth, are given in our illustrated Q-manual (Barton and 


Grimstad, 2014). Accelerated (i.e. increased) SRF values in the case of high stress give 


lower Q-values, and therefore the necessary heavier support (closer B c/c, thicker Sfr). 
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Competent rock, rock stress problems 
Low stress, near surface, open joints. > 200 < 0.01 


Medium stress, favourable stress condition. 
High stress, very tight structure. Usually 
favourable to stability, may be unfavourable for 


wall stability. 
Moderate slabbing after > 1 hour in massive rock. | 5-3 0.5-0.65 
Slabbing and rock burst after a few minutes in 0.651 
massive rock. 
Heavy rock burst (strain-burst) and immediate 
dynamic deformations in massive rock. 


Table 1. The 'accelerated' SRF value used in the Q-system S(fr) update of Grimstad and Barton, 1993 
when the estimated ratio of o,/0, reaches and exceeds 0.4. 


The common (but independently suggested) assumption by Grimstad and Barton, 
1993 and by Martin et al. 1998 that 0,/0,>0.4 will result in “stress-induced” fracturing, and 
the need for heavier tunnel support has recently been revisited. Research by co-author Shen 
using the fracture mechanics code FRACOD (Shen et al. 2013) shows that the assumed 
“high stress” fracturing is actually initiated by extensional strain, causing the weakest link 
(tensile strength 0,) to be exceeded first. Propagation by shearing may follow immediately if 
stress levels (or depth) are sufficient. Shearing may dominate at high stress or when a tunnel 
(or mine opening) is very deep, and even cause rock bursts, as propagation of fracturing in 
shear is unstable. (Shen and Barton, 2018). Figure 2 illustrates modes of fracturing seen 
or modelled in simulations of deep tunnels. Numerical modelling with the displacement 
discontinuity (DDM) based FRACOD code is also shown. This was a study specifically to 
see the effect of jointing on reducing the risk of rock-burst, due to dissipation of shear 
and tensile stresses with the help of the jointing. In deep rock tunnels, the presence of 
significant jointing appears to reduce the risk of fracturing and rock bursts, while in shallow 
tunnels, the presence of jointing, and of course faulting, increases the risk of collapse, if 
insufficient tunnel support is applied. The FRACOD models shown in Figure 2 show 1,000m 
deep simulations with O, «orizonta) = 2 Oo vertcay; The effect of jointing is to partly dissipate the 
fracturing of intact rock. (See numerous examples in Shen and Barton, 2018). Red fractures 
are tensile, green are shear, with respectively mode | and mode Il fracture toughness. In the 
top of Figure 2 we see a physical model of intersecting log-spiral shear fracturing caused by 
boring into a highly-stressed poly-axial (3D) flat-jack loaded cell of 50 x 50x 50cm size, in a 
direction deviated to one principal stress. (Addis et al., 1990). The three-dimensional nature 
of the fracturing, with fracturing also ahead of the face, has been experienced in some 
deep TBM tunnels. In the middle of Figure 2 (left) the world's first significant TBM tunnel is 
shown. This 2.2m diameter and nearly 1 km long tunnel was driven by a Beaumont TBM in 
1880, using steam-power. The chalk marl of the adjacent Channel Tunnels, driven between 
the UK and France some 110 years later, has a uniaxial compressive strength o, of only 4 
to 9MPa. The failure of the haunches in this pioneer pilot tunnel occurred where it passed 
under a 70m high cliff with a consequent approximate 2MPa increase of vertical stress. We 
can assume that the ratio of o,/0 ratio would have exceeded 0.4 by a significant margin. 
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Figure 2. Examples of 'stress-induced", or more correctly extension-strain induced fracturing, followed 
by propagation-in-shear. (Top) sandstone tunnel models, (middle) two real TBM tunnels from 1880 and 
2009 (chalk marl and marble). The two FRACOD models simulate deep TBM excavations: intact or 
jointed granite. 


The TBM photograph showing complete tunnel collapse is the result of a tragic rock- 
burst accident in a very deep TBM pilot tunnel of 5m diameter in China. This was bored in 
insufficiently strong, originally massive marble with o, of 70-120MPa. The tunnel depth was 
typically 1 to 2.5 km. All tunnels in this big hydropower project (Jinping Il) finally saw the 
TBM-driven headrace tunnels replaced by slightly less hazardous drill-and-blasting. In the 
latter, the highest tangential stresses are displaced a bit deeper into the surrounding rock 
due to the fracturing caused by the blasting. 
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FAILURE IN EXTENSION 


Based on the extension-strain theory, which was promoted by Stacey, 1981, if the 
strain in a given direction becomes tensile and reaches a critical value, tensile fracturing will 
occur. A two-dimensional equation for expressing extension strain (in the lateral direction) 


is as follows: 
e=lo,-vo]/E' (1) 


where v is the Poisson's ratio of the intact rock and E” is the generalised term for 
Young's modulus (E). (While E' = E for the plane stress condition; E' = E/(1-v?) for plane 
strain, i.e. when no expansion in third dimension. 

According to Equation (1), tensile strain may occur in a stress field where both 
principal stresses are compressive due to Poisson's effect. This explains why tensile 
fracturing can occur in the roof/wall of an underground opening, and also behind a cliff or 
mountain wall, where no tensile stress is expected. The only requirement will be that vo, > 
O, i.e. the disparity between the major principal stress (0,) and the minor principal stress 
(0,) needs to be high enough. This condition is met near to a tunnel wall, because the radial 
(= minor principal) stress becomes much lower than the tangential (= major principal) stress. 
Since Poisson's ratio is frequently about 0.25, a location where the stress ratio 0,,0,> 4 
is required as an absolute minimum, but in fact sufficiently diverging principal stresses are 
required to generate significant lateral strain to cause tensile fracturing. The critical tensile 
strain for tensile fracturing to occur can be determined using the tensile strength of the rock. 
We can simply express this as: 


e,=0/E' (2) 

Using the critical tensile strain in Equation 2 to replace e, in Equation 1 we obtain 
(by elimination of E” from both sides of the new equation), the critical compressive (i.e. 
tangential) stress for tensile fracturing (or spalling) to occur: 

Since 0/E'=[0,- vo ]/E' therefore: 

o,(spalling) = (0,+ 0,)/v (3) 

Considering that the confining stress o, is zero atthe wall of an underground opening 
(and next to a cliff or mountain wall) then for rocks with, typically UCS = 100, and Poisson's 


ratio = 0.25, it means that tensile fracturing will start when the tangential stress reaches = 
0.4 0,. This in fact is the simple arithmetic source of the “magic ratio” (0.4+0.1) that we saw 
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in Figure 1 following Martin et al. 1998, and confirming earlier observations of spalling and 
rock bursts in deep (or anisotropically stressed) Norwegian tunnels, as utilized by Grimstad 
and Barton, 1993 for quantifying SRF; the stress reduction factor, which is used to increase 
tunnel support via the rock mass classification Q-system. 


NEW FORMULAE FOR VERTICAL CLIFFS 


We will now leave the discussion of failure in deep tunnels, and concentrate on 
cliffs and mountain walls. We will also view evidence of extension strain fracturing in the 
form of planar sheet jointing from the world of rock climbing. Curved mountain slopes like 
the back of Yosemite's Half Dome in California are not the only reason for sheeting joint 
development. Figure 3 shows some compelling evidence for the range of “cliff' heights 
apparently caused by the wide range of compression strength of rock. The world's highest 
almost vertical mountain walls top out on either side of 1,300m. The probable record is 
1,340m for Great Trango Tower in the Karakoram, Pakistan. We may expect (laboratory- 
scale) compressive strengths of 100-150MPa for the granites in the highest walls, while in 
the case of Cappadocia's ancient cliff dwellings in porous tuff in Turkey, the strength may 
be only 1MPa, limiting cliff heights to the frequently observed 15-20m. The key to the huge 
range of heights illustrated in Figure 3 is that the tensile strengths of rocks may range 
from no more than 5-10MPa for exposed granite, down to 0.05-0.1MPa for exposed tuff. 
Numerous examples of the latter are given by Aydan and Ulusay, 2008. All tensile strengths 
tend to get downgraded by weathering during millennia of thermal cycles, which may exceed 
a 70-80ºC annual range, and vary strongly even during one day. The big range of what are 
in fact “limited' cliff and mountain-wall heights, can be roughly explained by comparing the 
assumed maximum vertical stress (the critical spalling stress in this case) with extensional 
strength 0,/v. To do this involves comparing cliff height and rock density, with the ratio of 
tensile strength and Poisson's ratio. 
= 100.0,/yv (4) 


critical 


(where q, is tensile strength in MPa, y is density in tons/mº, and v is Poisson's ratio). 
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Figure 3. The left-side photographs are examples of extreme mountain walls in hard or very hard rock. 
a) Great Trango Tower, Karakoram, Pakistan: approx. 1,350m, b) Mirror Wall, Baffin Island, Canada: 
approx. 1,200m, c) El Capitan granites, Yosemite: 950-1,000m. d) West Temple sandstones, Zion, Utah: 
650-700m, e) Beachy Head bedded chalk, England: 75-100m, f) Cappadocia tuff, Turkey: 15 to 20m. 
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Note that the multiplier of 100 in equation 4 is purely a function of convenience when 
using typical rock mechanics strength units, expressed in MPa. With strength in kN/m? and 
density in kN/mº the “100' can be dispensed with. With equation 4 in mind, one should refer 
to the worked examples in Table 2, where realistic values are employed, for a wide spectrum 
of possible rock conditions, ranging from uniaxial compression strengths of 150MPa typical 
for hard granite, down to just 1.0MPa which might be typical for a weakened porous tuff, as 
found in the cliffs of Cappadocia, depicted at the bottom of Figure 3. 


9c= 150MPa Sc = 50MPa Sc = 15MPa Sc = 1.0MPa 
Very hard Medium hard | Low strength Very low 
rock rock rock strength rock 


ao = 9MPa o = 3MPa o = IMPa o = 0.05MPa 
v=0.25 v=0.20 v=0.15 v=0.15 
o1/v=36MPa | o/v=15MPa | o/v=6.7MPa ot v = 0.33MPa 


y=28tonm” |y=2.5ton'm” |y=2.25tonm” |y=2.0ton'mº 
Hent= 1,290m | Henr= 600m Henit = 296m Horit= 16m 


Table 2. Examples of vertical height limits for near-vertical mountain walls and cliffs, over two orders of 
magnitude of rock strength, based on the application of equation 4. 


In general, the ratio of compressive and tensile strengths for rock are in the range 
of 5 < 0/0, < 20. However, in Table 2, since (compared to tunnels), cliffs and mountains 
are exposed to weathering forever, we made relatively conservative estimates of tensile 
strength o,. The proposed o/yv mechanism of extension failure has the effect of steepening 
and degrading mountain faces. The sketches and photographs presented in Figure 4 
show several aspects of the extensional mechanism. Because tensile strength is slowly 
reduced nearestto an exposed cliff or mountain face by constant cycling of temperature and 
moisture, it becomes easier for slabbing to occur. With assumed microcracking and grain- 
scale ice pressures, one may surmise an effective increase in Poisson's ratio, together with 
reduced near-surface tensile strength. If broadly correct, extension failure due to the 0,/v 
'mechanism” would become easier with time. This is probably the true cause of the relative 
frequency of slabbing due to the 'constant' propagation of sheet jointing. We will return to 
this subject later. 
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Figure 4. Top: Sketches illustrating the extensional strain failure mechanism involving the overcoming of 
the extensional strength 9, /v. Also shown is the possibility of rock avalanche-scale basal shear failures, 
due to adverse tectonic structures. This is slightly developed at El Capitan in Yosemite, 1,000m in height 
(bottom-left), but apparently well developed at Holtanna, a 750m high monolith, in Dronning Maud's 
Land, Antarctica (bottom-right). (Here there may be an unintended but probable camera rotation of 5º - 
the Holtanna potential failure plane” should probably be approximately 5º steeper.) 


The estimation of basal shear strength could be based on a combination of the 
strength components shown in Figure 5, since a final 'cascading' of failure may be involved 
(i.e. not linear Mohr-Coulomb or non-linear Hoek-Brown “c + o tan q”, but more correctly 'c 
then o tan q”, as increasingly emphasized in some recent literature). Most likely: failure of 
remaining intact 'bridges”, then shearing (or not) on these new, fresh, rough surfaces, then 
mobilization (or not) along the already established shear plane (or prior tectonic structure), 
and finally the limited strength of any clay-filled discontinuities or faults, displaying the 
lowest shear stiffness. (Barton, 2013). 
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or (308) 4, 


ra] 


(from Q-system) 


Figure 5. Representation of shear strength components in the form of laboratory tests, from Barton, 
1999. Several components of strength may be involved when a large-scale body of rock is approaching 
failure. These linear approximations are often misleading when a significant range of stress is involved, 

as emphasized by Barton, 2013, 2016. 


Although rock-surface curvature obviously helps for generating tensile stress as 
suggested by Martel, 2017, it is not a necessary condition for generating sheeting joints. 
The classic curved sheeting joints on 'the back of Yosemite's Half Dome have alternative 
means of development than the curvature, and indeed are remarkably planar up the 
750m vertical face. The planarity needs an explanation because sheeting joints are often 
completely planar over long distances, bot horizontally and vertically, as can be vouched for 
by rock climbers, and as can be seen in thousands of rock-climbing photographs. Interesting 
sources here are the following: Davis, 2013, Honnold and Roberts, 2016 and Florine and 
Moye, 2016. For the case of the Cappadocia tuffs, careful studies by Aydan and Ulusay, 
2003 and many others, have shown that temperature (freeze-thaw) cycling and moisture- 
content cycling have a degrading effect on compression and tensile strength. There is an 
effective annual temperature range of 70 to 80º C even in the shade, and this may even 
be increased for cliffs in direct sunlight. In the Alps, specifically the Matterhorn, the gradual 
degradation of the highest rock which is in almost permafrost conditions, has been nicely 
described by Weber et al., 2016. 


SOIL MECHANICS FORMULATIONS ARE NOT APPROPRIATE FOR FAILURE 
OF INTACT BRITTLE ROCK 

Concerning the question of limiting cliff and mountain wall heights it is appropriate 
to present and then check classic soil mechanics based solutions which at their simplest, 
involve a linear shear strength envelope defined by cohesion and friction (c and q). The 
Coulomb shear strength assumption, with allowance for Terzaghi's law of effective stress, 
assuming water pressure (u) along the failure plane, is as follows: 
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T=c+(0, -u)tany (5) 


Assuming a planar shear failure surface, dipping from the ground surface down to 
the toe of an imaginary vertical cliff, /owerbound and upperbound solutions can be obtained, 
based on limit equilibrium or limit theory analysis. Soils mechanics texts such as Verruijt, 
2001 indicate the following range of solutions for critical heights of vertical cuts in soils 
i.e.Fig. 6. 


2c/y tan(45º+0/2) < Hs 4c/y tan(45º+0/2) (6) 


Figure 6. Several approaches to the stability of a vertical cut appear in soil mechanics literature. The 
top example (a) shows the assumed equilibrium of three zones and gives a /owerbound solution. The 
bottom example (b) illustrates an upperbound solution involving a specific shear surface. After Verruijt 

(2001). The circular failure surface assumption was used by Fellenius, 1927. 
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Equation 6 gives a surprisingly large (2:1) range. Furthermore, because soil fails 
differently to intact rock, we find that the estimates of H, if (erroneously) used for rock cliffs 
or mountain walls are, remarkably, ffom 3 to 6 times in error, as indicated in the worked 
examples below. Note that a circular failure surface gives 3.85 for the multiplier. An exact 
solution is elusive according to Verruijt, 2001, and matters do not get easier if non-linear 
shear strength is considered. The lower bound of equation 6 is usually attributed to Drucker 
and Prager, 1952 and is referred to as the static solution, in contrast to the “dynamic” solution 
with multiplier “4”. 

Concerning the trial evaluation of “soil mechanics' formulations summarized in 
equation 6, we need to make appropriate estimates of cohesion. A lowerbound estimate 
of cohesion (c) for rock would be obtained most simply by assuming a straight-line, rather 
than a curved tangent between the uniaxial tension (0,) and uniaxial compression (0,) Mohr 
circles. The simple equation for the lowerbound cohesion intercept (c), derived from Mohr 
circle geometry, was given in Barton, 1976: 


c=Y(o.0)'? (7) 


We can examine the foregoing formulations with examples of a moderately strong 
valley-wall sandstone (UCS = 75 MPa) and a massive-scale mountain-wall in granite (with 
UCS = 150MPa). We will assume that 0,/0,= 15. (The expected laboratory-test range of 0,/ 
o, could range from 5 to 20: but 8 to 16 is more likely). The following strength assumptions 
are used: 


1. Sandstone o = 75MPa, o, = 5MPa c = Y(75x5) '2=9.7MPa 
2. Graniteo = 150MPa, o,= 1OMPa c = Y(150x10) "2=19.4MPa 
The gradient (p = internal friction angle) of the presently assumed straight line 


between the tensile and compressive strength Mohr circles, giving a /owerbound value of c 
is as follows: 


o/a =tan?(45º+0/2) (8) 


A friction angle & of approximately 61º is indicated if 0,/0, =15, as assumed. Using 
justthe lowerbound soil mechanics solution given by equation 6, serious errors of mountain 
wall heights are evident. Note the following density assumptions: 25kN/mº for sandstone, 
27.5kN/mº for granite. 

Sandstone valley wall" H,= 2c/y tan (45º+4/2) H =2x9.7x1000/25 tan (45º + 30.5º) 
= 3,001m 


Granite 'mountain-wall"H =2x19.4x1000/27.5 tan (45º+ 30.5º) = 5,456m 
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By comparison equation 4 (specifically applying to brittle rock) suggests H limits of: 
Sandstone: 100.0,/yv =100.5/2.5 x 0.25 =800m 


Granite: 100.0,/v =100.10/2.75 x 0.25=1,456m 


We may note from these solutions for rock, using the new formulation involving 
extensional theory (equation 4) that the real cases depicted in Figure 3 are informing us 
that it is not realistic to reckon with “laboratory scale” (i.e. unweathered, optimal samples) 
tensile strengths. These two mountain-wall estimates of 800 and 1,456m are perhaps 20- 
30% too high, but much better than 300 or 600% too high if erroneously trying to apply soil 
mechanics methods (equation 6). Note again the more sensible range of height limits given 
in Table 2. 


PLANAR SHEETING JOINTS FROM EXTENSION STRAIN MECHANISMS 


To conclude this section about a new method of estimating cliff and mountain-wall 
heights, we can tentatively apply the o, /v method to explain the origin of sheet jointing. 
Our seemingly almost fearless free-solo rock climbers, who are constantly climbing these 
planes, and finger-wedging up the related sub-vertical crack systems. Statistics from well- 
documented climbing routes in the Yosemite Valley given by Stock et al. 2012 indicate that 
'slabbing' or the fall of loosened sheeting joints is quite frequent. Indeed, there are frequent 
reports of changed climbing routes due to their occurrence, perhaps as frequently as every 
decade. 

Figure 7 is a good demonstration of the extreme planarity and smoothness that can 
be a common characteristic of sheet-jointing or exfoliation, in this case in Zion sandstone 
and the vertical (front) wall of Yosemite's (Half Dome) granite. A curved surface, facilitating 
the tensile component, is not in fact a necessary condition for the development of sheet 
jointing, as argued in a comprehensive (two centuries) review by Martel, 2017. Note that 
the long sub-vertical cracks loved and indeed needed by rock climbers, are probably a 
large-scale expression of extensional fracturing in the perpendicular direction. An excellent 
example is shown in Figure 8, with Alex Honnold again free-soloing. 
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Figure 7. Two famous free-solo rock climbers: Steph Davis and Alex Honnold on sheeting joints in Zion 
national park and in Yosemite national park (“Thank God Ledge, Half-Dome) where Honnold made the 
first ever free-solo ascent. Refer to Davis, 2013, and Honnold and Roberts, 2016 for some remarkable 
examples of planar sheeting joints. Rock and mountain climbing, as shown in countless internet sites, 
is a particularly rich source of examples of rock exposures at all scales, for those interested in rock 
fracture mechanics principles. 


MAXIMUM SLOPE ANGLES WHEN ROCK MASSES ARE JOINTED 


This paper has so far treated rock masses as if they were intact or sparsely jointed, 
both in the case of the fracturing in deep tunnels, and in the case of the limiting heights of 
cliffs in weaker rock, and mountain walls in strong rock. If instead, jointing is present, it has 
been assumed noto significantly alter the extension strain fracturing mechanism. However, 
is doubtful that this assumption is generally valid, but it may be acceptable in the case of 
cliff failures in (horizontally) bedded rock, since the tensile failure planes would be at right- 
angles to such beds. See for instance the cliff-front failures in chalk, Figure 3e. 
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Figure 8. Free-solo climber Alex Honnold on El Capitan. Extension fractures in the third dimension are 
the assumed origin. Such features may be hundreds of meters in extent and seem to be formed by 
mountain-induced stress (and strain) rather than being a part of a pre-existing major joint pattern. 


Methods have recently been developed for selecting safe and maintenance-free rock 


slope angles in more general jointed conditions. The method is called Q with general 


slope” 
similarity to the Q-system for tunnels. In fact the first four parameters are Enade except 
for the use of oriented Jr/Ja ratios to allow for wedges formed with two joint sets having 
different Jr/Ja magnitudes, and therefore potentially different influences on instability. The 
method was introduced by Barton and Bar, 2015 and more recently described by Bar and 
Barton, 2017 after the collection of more than 400 case records, mostly by co-author Neil 
Bar. 
Oslope = a 2r e ouico (9) 
dy Ja 0 SR ope 


As with the Q-system, the rock mass quality in Q-slope can be considered a function 
of three parameters, which are crude measures of: 
1. Block size: (RQD/ ). 


2. Shear strength: (J /J) or average shear strength in the case of wedges (J /J ), 
X (S/S) 
3. External factors and stress: (J,.,/SRF supe): 
Barton and Bar, 2015 derived a simple formula for the steepest slope angle (B) not 
requiring reinforcement or support for slope heights less than 30m. This formula has now 
been extended to much larger slope heights: 
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B = 2010819 Oyjope + 65º (10) 


Equation 9 matches the central data for stable slope angles greater than 35º and 
less than 85º. From the Q-slope data, the following correlations are simple and easy to 
remember: 


* - Q-slope = 10 - slope angle 85º 
* Q-slope=1 -slope angle 65º 
* - Q-slope=0.1 - slope angle 45º 


*- Q-slope = 0.01 - slope angle 25º 


Numerous case records are illustrated in Figure 9. The following example from Bar 
and Barton, 2017 gives some insight into this simple method, and Figure 10 illustrates the 
concept for the case of open pit bench angles, obviously in jointed, as opposed to sparsely 
jointed rock. 

A 30m high slope was excavated at an angle of 65º and failed shortly after. The 
wedge failure occurred in weak, moderately weathered sandsione (o, = 35MPa). The 
following Q-slope ratings were assigned during the back-analysis: 


RQD = 40-50% 
J,=9 
SetA:J =1,J,=4, O-factor = 0.5 


Quioye Stabllity Chart 


s 


g 


— 0Q-Slope Equaton 


8 


Slope Angle, f (degrees) 
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Quupo Case Study Data 
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Figure 9. The coloured areas indicate stable (green) and unstable (red). The case records show fairly 
consistent separation, with a transition zone (grey). Bar and Barton, 2017. 


SetB:J,=3,J, = 4, O-factor = 0.9 
Set C: Release plane or tension crack that did not contribute to the overall shear 
strength of the wedge. 


J...= 1 (desert environment, competent rock and generally stable structure where 
Set B has limited continuity). 


SRF, = 2.5 (slight loosening due to surface location), SRF, = 2.5, SRF = N/A. 
Based on the assigned ratings, Q-slope and B were estimated as follows: 


ssi(1 3 l 
=""4| -x0.5 | — x0.9 | |x — = 0.206 
Osiope 9 da Vá: Je Zs 


B = 2010g,9(0.206) + 65º = 51º 


wice 


=0.14 
Bench Face Angle = 45º 


Quipe=.42 
Bench Face Angle = 55º 


Quioçe=0.94 
Bench Face Angle = 65º 


Figure 10. An open-cast slope in Laos, showing increasing values of Q,,, and slope angles, as greater 
depth and sounder rock is reached. Bar and Barton, 2017. 


Q-slope suggested an angle of 51º would have resulted in a stable slope (i.e. 
approximately 15º shallower than excavated and consistent with kinematic analysis). 
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MOUNTAIN HEIGHTS LIMIT OF 8 TO 9KM 


Increasing the scale dramatically, one can tentatively suggest that the highest 
mountains of 8,000 and 9,000m are not due to the “limited” uniaxial compression strength as 
sometimes proposed, but are due to critical state, non-linear rock mechanics. The limiting 
strength is more likely to be given by the top (horizontal) part of a rock's strongly curved 
shear strength envelope, or perhaps by the slightly lower brittle-ductile transition. Figure 11 
in fact demonstrates similar magnitude for the critical confining pressure and the magnitude 
of uniaxial compressive strength, such as 200 MPa, as shown by Singh et al. 2011. 

Note that the overall curvature of the suggested shear strength envelope (Figure 
11) is somewhat greater than that of Hoek-Brown, which is presently one of the most used 
non-linear strength envelopes for intact rock. 

The world's highest mountains (14 peaks of 8,000 to almost 9,000m height, with 
Everest at 8,848m: Figure 12) cannot be limited by the uniaxial compression strength of 
rock. This is because the rock strength would need to be the confined compression strength 
(o, in Figure 11), since great depth is obviously involved where stress and strength are 
(almost) in equilibrium. The possible range of perhaps 600 to 900MPa for the confined 
compressive strength of strong igneous or metamorphic rocks at 10km depth would support 
mountains of 20 to 30 km height. 


3 
7 
: 
o 
=== TRANSITION 
[o] 
E 
V E 
> 


x 
1= Uniaxial Tension >= - 
2 = Uniaxial Compression 

3 = Brittle-Ductile Transition 


4= Critical State 3 353 


Figure 11. The curved failure envelope for intact rock up to the critical state proposed by Barton, 1976. 
The critical state marks the location of the maximum possible shear strength for a given rock type. 
The uniaxial tensile (1) and uniaxial compression (2) Mohr circles are shown on the left, and provide 
the lower-bound estimate of cohesion c, that was utilized earlier. Singh et al. 2011, in a major review, 
have shown that the critical confining pressure is equal or close to, the uniaxial compression strength 
(see small blue arrow). Both the deviation from linear Mohr-Coulomb, and the equation of the curved 
envelope are developed (by Singh et al. 2011, and by Shen et al. 2018). 
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The empirical evidence of millions of years is therefore violated if trying to use the 
(confined) compression strength of the rock, and rock mechanics principles are actually 
violated by those suggesting use of the 200-300 MPa if referring to the use of uniaxial 
strength. The fourteen mountains in the exclusive 8 to 9 km height class provide empirical 
evidence of total vertical stresses that might be as high as 9,000 x 3/100MPa = 270MPa, 
assuming a composite density as high as 3.0. If there could be a pore/joint water pressure 
as high as 5OMPa (with some upper drainage into valleys assumed here), then tentative 
application of the law of effective stress brings us to a principal effective vertical stress of 
220MPa. This would apply in the lower part of a potential shear failure surface. However, we 
need to consider the shear stress generated by this effective vertical stress of 220MPa. On 
a nominal plane inclined at 45º this could be as low as 150MPa, considering round figures. 
In case of an elevated horizontal stress (i.e. k, > 1), equilibrium would be improved, and a 
higher mountain could be supported before reaching a shear strength limit. 


Figure 12. Mount Everest at 8,848m (extract from Wikipedia photo). It is proposed that it's maximum 
height is shear-strength limited, following a critical state rock mechanics limit, as shown in Figure 8. 


CONCLUSIONS 


Tensile strength and Poisson's ratio explain the limited maximum heights of cliffs 
and steep mountain walls, and the origin of planar sheeting joints. A range of maximum 
heights from 20m in tuff, 100m in chalk, 650m in sandstone, to 1,300m in granite can be 
sensibly quantified by considering failure caused by extensional strain and fracturing in 
tension in each case. There are parallels in the world of deep tunnels in hard rock. The 
widely quoted critical tangential stress of 0.4 x UCS that may be reached by deep hard-rock 
tunneling should be replaced by the ratio o,v, i.e. initial tensile failure which is mobilized 
by extensional strain. These two ratios are numerically equivalent. Shear strength and 
tensile strength, ably assisted by Poisson's ratio, are inevitably the weakest links in “high- 
stress' structural geology and in the more typically low-stress' processes in geomorphology, 
respectively. When rock masses are jointed the extension strain-induced failure mechanism 
may no longer apply, and stability has to be assessed by characterizing the properties of 
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the jointing. The Ou methodis proving to be useful here. The highest mountains of 8,000- 
9,000m are limited by maximum shear strength, not by compressive strength. The confined 


strength of competent mountain-forming rock is several times too high. 
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RESUMO: A geração da drenagem ácida de 
mina (DAM) ocorre a partir da oxidação de 
sulfetos em rochas mineralizadas, sendo a 
pirita o sulfeto mais comum nestas rochas. 
ADAM pode ter pH's variando de 1a 6e em 
virtude disto possui um elevado potencial de 
dissolução de minerais e de solubilização 
e transporte de elementos tóxicos (As, 
Cd, Cu, Ni, Pb, Cr, Hg, Zn) presentes 
rochas mineralizadas, contaminado solos, 
rios e aquiferos. No Brasil, as principais 
ocorrências de DAM descritas são a mina 
de urânio em Caldas, de cobre em Carajás 
e nas minas de carvão na região sul. A 
predição de geração e neutralização da DAM 
é feita por experimentos com os rejeitos de 
minérios e rochas estéreis associadas por 
meio de experimentos cinéticos e estáticos. 
A modelagem geoquímica é uma ferramenta 
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com potencial de utilização para a avaliação 
teórica da capacidade de geração e 
neutralização da DAM antes mesmo 
da realização de experimentos. Como 
estudo de caso foram realizados ensaios 
laboratoriais de geração e neutralização da 
DAM em estéreis da mina de urânio Osamu 
Utsumi: (i) potencial máximo de acidez, (ii) 
teste geração de acidez liquido de adição 
única e (iii) avaliação da neutralização 
da ácidez - BC Research. Foi realizada a 
modelagem geoquímica com o software 
PHREEQC de todos os experimentos. 
Por meio da modelagem geoquímica foi 
possível reproduzir os ensaios laboratoriais 
e predizer as capacidades de geração e 
neutralização da acidez. Além disto, por 
meio da modelagem geoquímica foi possível 
estimar a quantidade solubilizada de 
minerais, responsável pela baixa capacidade 
de neutralização das rochas. Desta forma, a 
modelagem geoquímica se mostrou eficaz 
na predição da geração e neutralização da 
DAM, podendo ser utilizada como ferramenta 
exploratória das condições ambientais antes 
mesmo da realização dos experimentos, 
a partir do conhecimento detalhado da 
mineralogia das rochas. 

PALAVRAS-CHAVE: Drenagem ácida de 
mina, geração, neutralização, modelagem 
geoquimica, PHREEQC. 
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GEOCHEMICAL MODELING AND EXPERIMENTAL TESTS TO DETERMINE 
THE GENERATION AND NEUTRALIZATION POTENTIALS OF ACID MINE 
DRAINAGE: CASE STUDY OF THE OSAMU UTSUMI URANIUM MINE (CALDAS- 
MG) 

ABSTRACT: The generation of acid mine drainage (AMD) occurs from the oxidation of 
sulfides in mineralized rocks, with pyrite being the most common sulfide in these rocks. AMD 
can have pH's ranging from 1 to 6 and as a resultit has a high potential for dissolving minerals 
and solubilization and transport of toxic elements (As, Cd, Cu, Ni, Pb, Cr, Hg, Zn) present. 
mineralized rocks, contaminated soils, rivers and aquifers. In Brazil, the main occurrences 
of AMD described are the uranium mine in Caldas, the copper mine in Carajás and the 
coal mines in the southem region. The prediction of generation and neutralization of amd is 
made by experiments with ore tailings and associated waste rocks through kinetic and static 
experiments. Geochemical modeling is a tool with potential use for the theoretical evaluation 
of the generation and neutralization capacity of AMD even before carrying out experiments. As 
a case study, laboratory tests were carried out on the generation and neutralization of AMD in 
waste from the Osamu Utsumi uranium mine: (i) maximum acidity potential, (ii) single-addition 
liquid acidity generation test and (iii) evaluation of the neutralization of acidity - BC Research. 
Geochemical modeling was carried out with the PHREEQC software for all experiments. 
Through geochemical modeling, it was possible to reproduce the laboratory tests and predict 
the acidity generation and neutralization capabilities. Furthermore, through geochemical 
modeling it was possible to estimate the solubilized quantity of minerals, responsible for 
the low neutralization capacity of the rocks. In this way, geochemical modeling proved to 
be effective in predicting the generation and neutralization of AMD, and can be used as an 
exploratory tool for environmental conditions even before carrying out the experiments, based 

on detailed knowledge of the mineralogy of the rocks. 
KEYWORDS: Acid mine drainage, genaration, neutralization, geochemical modeling, 
PHREEQC. 


INTRODUÇÃO 


A predição de drenagem ácida de mina (DAM) é um importante método para 
prevenir ou mitigar possíveis impactos ambientais causados pela oxidação de minerais 
sulfetados. A DAM ocorre a partir da oxidação dos sulfetos presentes na pilha de rejeitos, 
gerando ácido sulfúrico, reduzindo o pH da solução para valores entre 1 e 6 (INAP, 2009). 
Com a diminuição do pH do meio, elementos tóxicos como As, Cd, Cu, Ni, Pb, Cr, Hg, Zn, 
são dissolvidos e tem sua mobilidade magnificada, contaminando corpos hídricos, solos 
e podem ser incorporados na cadeia alimentar de plantas e animas, incluindo o homem 
(Figueiredo, 2000, Blowes et al.,2014). Lottermoser (2007) descreve a reação (1) de 
formação da DAM a partir da pirita (Fes,). 


Fes, (s) + : 0, (9) + H,0 (1) > Fe?* (aq) + 250,” (aq) + 2H* (aq) (1) 
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O processo de oxidação de materiais sulfetados e a formação de DAM envolvem 
processos químicos e biológicos e podem variar a depender das condições climáticas, 
geológicas e ambientais (INAP, 2009, Nordstom & Alpers, 1999). A influência de bactérias 
do gênero Thiobacillus aumenta a taxa da reação significantemente em relação a taxa de 
reação sem a presença de microrganismo (Robertson & Broughton, 1992). 

No Brasil, as principais ocorrências de drenagem ácida estão na mineração de 
urânio em Caldas — MG, cobre em Carajás — PA e nas minas de carvão na região sul do 
país (Borma, 2006). Em minas de urânio além da própria drenagem ácida há o agravante 
da presença dos elementos radioativos, em virtude disto, países como Austrália, Canadá, 
Estados Unidos e França têm desenvolvido exaustivos estudos de impactos da drenagem 
ácida de minas de urânio em operação e desativadas (Cipriani, 2002). 


Métodos para avaliação de drenagem ácida de mina 


Os métodos de avaliação e predição da drenagem ácida de mina são constituídos 
por experimentos de rejeitos de minérios e rochas estéreis associadas com o intuito de 
se avaliar o potencial de gerar ácido, o potencial de neutralização destas rochas. Estes 
métodos são denominados cinéticos e estáticos. Os métodos cinéticos predizem o potencial 
de geração de drenagem ácida em longo prazo, podendo ser analisado pH, Eh, acidez, 
alcalinidade, condutividade elétrica, metais dissolvidos, taxa de oxidação e produção de 
sulfato (U.S.EPA 1994). Esses métodos se aproximam das condições climáticas e dos 
processos encontrados na localização da mina, com uma taxa acelerada. (Farfan, 2004). 
Os métodos estáticos se baseiam em avaliar o equilíbrio entre a geração de ácido e a 
capacidade de neutralização da acidez, caraterístico como equilíbrio ácido base, do inglês 
“Acid Base Accounting” (ABA). 

Arealização de diferentes ensaios para a avaliação de DAM permite uma comparação 
dos resultados alcançados, oferecendo uma maior confiabilidade dos dados obtidos, e 
a comparação dos resultados, custos e dificuldades de execução dos testes, ajuda na 
escolha de métodos mais objetivos para a predição, mitigação e remediação da DAM. 

Além dos métodos experimentais a modelagem geoquímica tem sido utilizada com 
sucesso na predição e remediação da geração de DAM, na simulação de migração de 
plumas contaminantes e na avaliação da atenuação natural. A modelagem geoquímica 
pode prover insights em condições futuras; determinar quais variáveis são as mais 
importantes na determinação das condições futuras; avaliar os efeitos das abordagens 
alternativas para o gerenciamento da DAM; avaliar efeitos potenciais de parâmetros 
incertos; estabelecer objetivos e condições de teste para estudos de campo e laboratório e 
integrar as informações disponíveis (Martin et al., 2005; Blowes et al., 2014; Nordstrom & 
Campbell, 2014; Nordstrom & Nicholson 2017; Verburg et al., 2009). 
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O objetivo deste trabalho é comparar os ensaios experimentais de laboratório acerca 
da capacidade de geração e de neutralização da DAM e a predição destes experimentos 
pela modelagem geoquímica, tendo como estudo de caso os rejeitos do Bota Fora 4 da 


mina de urânio Osamu Utsumi. 


A mina de Osamu Utsumi 


A mina de Osamu Utsumi se situa sobre as rochas ígneas da região de Poços 
de Caldas. A caldeira de Poços de Caldas tem forma quase circular, com diâmetro de 
aproximadamente 30 km. Associadas a estas rochas ocorrem mineralizações de urânio, 
molibdênio, tório e terras raras (Holmes et al. 1992 e Fraenkel et al. 1985). Situado numa 
região de transição entre a Bacia Sedimentar do Paraná e a Serra da Mantiqueira, o 
Complexo Alcalino de Poços de Caldas é constituído por rochas sedimentares referida 
bacia, embasadas por rochas magmáticas e metamórficas do Complexo Pinhal (Wernick & 
Penalva 1980). O Complexo Alcalino de Poços de Caldas resulta de uma atividade ígnea 
policíclica, que inclui processos metassomáticos intensos e um forte intemperismo. Estes 
processos originaram vários litotipos e mineralizações de interesse econômico, como de 
minério urano-molibdenífero, zircono uranífero, zirconífero e as bauxitas (Fraenkel et al. 
1985). As rochas ígneas alcalinas do Complexo Alcalino de Poços de Caldas pertencem 
ao grupo nefelina-sienitos e seus correspondentes afaníticos (tinguaítos e fonólitos). As 
rochas alcalinas hidrotermalizadas, encaixantes da mineralização uranífera, evidenciam 
empobrecimento em sílica e alumínio e um forte enriquecimento em potássio, sódio e cálcio. 
Destacam-se neste complexo alcalino rochas com altos teores de titânio, com minerais de 
molibdênio, zirconita, pirita e fluorita (Fraenkel et al.1985). Na Jazida do Cercado, local 
onde foi lavrada a mina Osamu Utsumi, o urânio é encontrado sob forma de uraninita 
[UO,] e coffinita [U(SiO4)1-x(OH)4x]. A uraninita se encontra normalmente em brechas 
que contêm zirconita (ZrSiO4), baddeleíta (ZrO02), jordisita (MoS2), pirita (FeS2), esfalerita 
(ZnS), galena (PbS) e fluorita (CaF2), sendo que argilas também podem estar associadas 
a uraninita maciça em fraturas (Oliveira 1974, in Tedeschi 2005). 


Passivo ambiental 


A mina se manteve em atividade entre os anos 1982 e 1995 para extração de urânio. 
Após o término das atividades de extração, dentre os inúmeros passivos ambientais de 
responsabilidade das Indústrias Nucleares do Brasil (INB), a geração da DAM a partir de 
pilhas de estéril tem sido o de maior impacto. Este passivo é o que tem demandado maiores 
preocupações e tratamentos em virtude da possibilidade de contaminação de drenagens 
e águas subterrâneas que extrapolem a área da empresa, alcançando a população 
circundante (Alberti, 2017). O material estéril retirado da cava da mina foi depositado em 
pilhas a céu aberto denominadas bota-fora (BF) localizados no seu entorno (BF1, BFS, 
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BF4, BF7 e BF8) e no seu interior (BF cava). O material estéril com característica terrosa, 
removido das camadas superiores da área lavrada foi depositado dos bota-fora BF1, BF2, 
BF3 e BF7. Os materiais com característica rochosa, extraídos nas camadas inferiores, 
foram depositados nos bota-fora BF4, BF8 e BF cava. O BF 4 (Fig. 1) possui uma área de 
56.4 ha e contém 12.4 mil m3 de estéril, e o BF8 tem uma área de 64.4 ha e 15.0 milhões 
m3 de estéril (Cipriane, 2002; Alberti, 2017). Os estéreis do BF4 e BF8 são constituídos 
por rochas tinguaíticas/ fenolíticas, onde minerais dominantes são o feldspato potássico 
e a sericita, e subordinadamente minerais associados ao hidrotermalismo como pirita, 
fluorita, minerais de urânio, molibdênio e zircônio (Cipriani 2002). A DAM gerada a partir da 
oxidação da pirita, é tratada por meio da adição de cal hidratada. Em virtude do tratamento 
há a precipitação de particulados em bacias de decantação. Esse material particulado é 
depositado na cava da mina e os efluentes são liberados no córrego do Cercado (Cipriani, 
2002). 

A DAM gerada no BF4 (Souza et al., 2016) possui pH 3.77 e Eh 457.9 mV, condições 
que elevam a concentração de metais na forma solúvel como o Al (122.3 mg L-1), U (4.56 
mg L-1), Th (0.099 mg L-1), além do metaloide As, na forma de H,AsO, (0.90 mg L-1). 
As águas oriundas da pilha de BF4 seguem para a Bacia Nestor Figueiredo, e então 
bombardeadas para a área de tratamento AA 440. (Freitas e Chapadeiro, 2012). A DAM 
gerada no BF8 apresenta pH 3.80 e C.E. 900 yS/cm, concentrações elevadas de Al (52.60 
mg/L), Mn (39.8 mg/L), Pb (0.09 mg/L), Zn (4.33 mg/L), Fe (1.73 mg/L) F (24.1 mg/L) e U 
(1.05 mg/L). A bacia Ivan Antunes recebe as águas infiltradas do BF8, então segue para o 
tratamento na área AA 440. (Freitas e Chapadeiro, 2012). 


a b 
Figura 1: Foto do BF4. A - Vista de cima do talude. B - Vista de baixo do talude. 
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MATERIAIS E MÉTODOS 


Amostragem 


Para a realização dos ensaios de geração e neutralização de acidez foram coletadas 
amostras de rochas e materiais sólidos em pelo menos 10 pontos do BF4, de modo a se 
ter uma representatividade e variabilidade espacial. Estas amostras foram caracterizadas 
mineralógicamente. Foram coletados cerca de 37 kg de blocos de rocha soltos no talude 
do BF4. As amostras, foram lacradas e identificadas em sacos plásticos e levadas para o 
laboratório. 

No laboratório as amostras foram separadas e classificadas perante suas 
características macroscópicas, e com base no NAGRA (Waber, 1990; Schorscher, 1991) da 
mina. As rochas foram classificadas em dois grupos (Fig. 2): o grupo brechoide, composto 
por rochas de granulação fina com coloração cinza clara, correspondendo a 34.1 kg de 
amostra. O outro é grupo é denominado de Ferruginoso, com 3.1 kg, e composto por uma 
rocha fina com uma carapaça de rígida de óxido de ferro. 


dl 


6 
5 
4 
3 
2 


Figura 2: a - Amostra do grupo “Ferruginoso”. a - Amostra do grupo “Brechoide”. 


Os grupos de rochas foram britados separadamente por um britador de mandíbulas e 
peneirados em diferentes tamanhos de malha. Foram obtidas as classes entre os intervalos 
do tamanho da malha, conforme a Tabela 1, em ABNT e milímetros. A classe “Total” é a 
junção das demais classes com composição proporcional a porcentagem obtida, conforme 
a Tabela 2. 


ABNT >4 4-10 10-16 16-35 35 - 60 < 60 Total 
mm >4,76 4.76 - 2.00 1.19 1.19-0.50 0.50- 0.25 <0,25 Total 
2.00 


Tabela 1: classes a partir do tamanho da malha da peneira. 


Descobertas e desafios nas geociências Capítulo 4 


60 


CLASSES BRECHOIDE FERRUGINOSA BRECHOIDE FERRGINOSA 


(KG) (KG) (%) (%) 

>4 7.2 0.6 21.2 18.4 
4-10 10.5 1.0 30.7 31.5 
10-16 3.3 0.3 9.7 10.3 
16-35 5.5 0.5 16.0 15.2 
35-60 2.8 0.2 8.1 77 
<60 4.9 0.5 14.3 16.9 
TOTAL 34.1 3.1 100 100 


Tabela 2: Massa e porcentagem das classes dos dois grupos de amostras 


A escolha do tamanho da malha das peneiras foi devida a preparação das amostras 
descrita em Method 1627: Kinetic Test Method for the Prediction of Mine Drainage Quality. 

Por último foi realizada a moagem das amostras classificadas no peneiramento em 
um moinho de disco de ágata no laboratório de preparação de amostras do IG/UNICAMP. 
Aproximadamente 150 g do material de cada classe dos dois grupos foi homogeneizado, 
quarteado e por fim moído. 


Caracterização mineralógica 


As rochas dos dois grupos foram moídas em moinho de ágata até a granulometria 
inferior a 0.063 mm. Posteriormente as amostras foram analisadas por difração de raios 
X no Laboratório de Difração de Raios-X do IG/UNICAMP. Além da determinação dos 
minerais foi realizado sua quantificação por meio da modelagem Rietveld. 


Avaliação da geração de drenagem ácida 


POTENCIAL MÁXIMO DE ACIDEZ (AMIRA, 2002) 


O potencial máximo de acidez (PMA) é definido pela quantidade de enxofre (S) na 
amostra. O cálculo assume que o enxofre quantificado ocorre como pirita (FeS,), onde 
reage em condições oxidantes gerando ácido, conforme a reação 2. 


15 7 
FeS, 0a +5H0 > Fe(0H), + 2 H,S0, (2) 


De acordo com a reação 2, 1% de S como pirita é equivalente a 30.6 kg de H,SO, 
por tonelada de material. O PMA pode ser calculado utilizado a equação (1) abaixo: 
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TESTE GERAÇÃO DE ACIDEZ LIQUIDO DE ADIÇÃO ÚNICA (AMIRA, 2002) 


O teste NAG (do inglês “Net Acid Generation” é realizado com a finalidade de se 
contabilizar o potencial líquido de geração ácida. Envolve a reação de uma amostra com 
peróxido de hidrogênio para oxidar rapidamente quaisquer minerais sulfetados contidos em 
uma amostra. As reações de geração de ácido e de neutralização de ácido podem ocorrer 
simultaneamente. Portanto, o resultado representa uma medida direta da quantidade 
líquida de ácido gerada pela amostra. 

São oxidados 2.5 g de amostra pulverizada (< 74 um) com uma solução de 15% 
peróxido de hidrogênio H,O, (30% de H,O, diluído 1:1 com água destilada). Deixa-se então, 
a reação ocorrer por aproximadamente 12 horas. A reação é então aquecida em uma placa 
quente até a efervescência parar ou num mínimo de duas horas. Após o resfriamento 
à temperatura ambiente é medido o pH da solução, nomeado de NAGpH, então as 
suspensões são tituladas até pH 4.5 com solução de NaOH 0.1 mol/L padronizada. O 
potencial de geração ácida é dado pela equação 2 abaixo. 


49 xV xM 
NAG = idas (2) 
m 
NAG = potencial de geração ácida (kg H,SO tonelada) 
V = volume de NaOH usado na titulação (ml) 
M = concentração de NaOH usado na titulação (moles/litro) 


m = massa da amostra reagida (9) 


O resultado do NAG é interpretado a partir da Tabela 3 abaixo. 


NAGpH NAG (kg H,SO,4) Potencial ácido da amostra 
=4.5 0 Não formação ácida 
<4.5 =5 Potencial de formação ácida baixa capacidade 
<4.5 Ex) Potencial de formação ácida 


Tabela 3: Interpretação dos resultados do teste NAG. 


Descobertas e desafios nas geociências Capítulo 4 


62 


Avaliação da neutralização da ácidez 


B.C RESEARCH (DUNCAN & WALDEN 1975) 


Esse método tem como objetivo mensurar a quantidade de bases presentes na 
amostra titulando uma suspensão da amostra com ácido sulfúrico a pH 3.5. Para isto, 
aproximadamente 10g de amostra pulverizada são suspensas em 100 ml de água deionizada. 
As suspensões, então, são agitadas por 15 minutos. Posteriormente as suspensões são 
tituladas com solução padrão de H,SO, 0.0103 mol/L até pH 3.5, tomando-se o cuidado 
para que este valor não se altere ao longo de 4 horas. O potencial de neutralização é dado 
pela equação 3. 


CxVxM 
= m 


PN 


PN = potencial de neutralização (kg H,SO tonelada) 
C = concentração de H,SO, (mol/L) 

V = volume de H,SO, (ml) 

M = massa molecular do H,SO, 

m = massa da amostra reagida (9) 


Modelagem geoquimica 

A modelagem geoquímica foi realizada utilizando-se o software PHREEQC 
(Parkhurst & Appelo, 2013) com o banco de dados lIlnl.dat. Os minerais utilizados foram 
K-feldspar (como microclínio), illite (como fengita) e a pirita. Foram simulados todos os 
procedimentos experimentais de determinação de geração e neutralização de DAM. 


RESULTADOS E DISCUSSÕES 


Caracterização mineralógica 


Os principais minerais identificados nas amostras por meio da difração de raios x 
foram microclínio, pirita e fengita. Para os experimentos e modelagem geoquímica, em 
virtude da homogeneidade mineralógica foi utilizada a classe Total do grupo Brechoide, por 
ser representativa de todo rejeito. Os teores dos minerais presentes na rocha, obtidos por 
modelagem Rietveld são apresentados na Tabela 4. 
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Y% em massa Mol Mols em 1 kg de rocha 


microclínio 82.4 278 218) 
fengita 14.1 389 0.3 
pirita 3,4 120 0,3 


Tabela 4: Composição mineralógica da rocha em % em massa obtida pela modelagem Rietveld. 


Potencial de geração ácida 


POTENCIAL MÁXIMO DE ACIDEZ 


O PMA é calculado a partir da quantidade de enxofre presente na amostra utilizando 
a equação (1). A quantidade de enxofre (%S), em porcentagem, é calculada a partir da 
proporção entre a porcentagem de pirita (3,4%), encontrada a partir da modelagem Rietveld 
no tópico acima, e a massa molecular da pirita (MMp) e 2 vezes a massa molecular do 
enxofre (MMS), a pirita (FeS,) contém 2 átomos de enxofre (S). A equação (4) representa o 
cálculo feito para calcular a porcentagem de enxofre na amostra. 


34% x MMS 
W=————— 
ad MMp (4) 
MMS = 64,120 (g/mol) 


MMp = 119,975 (g/mol) 


A porcentagem de enxofre na amostra encontrada é de 1.814%. Portanto para a 
equação (1) é calculado um PMA de 55.51 kg de H,SO,/ tonelada de amostra. 

A modelagem geoquímica do Potencial Máximo de Acidez foi realizada levando em 
consideração a proporção dos minerais presentes em 1 kg de rocha em equilíbrio com água 
pura, num sistema aberto à atmosfera e sem permitir a precipitação de nenhum mineral. 
No equilíbrio químico houve a solubilização total da pirita e parcial da ilita, a concentração 
do microclínio não foi afetada. O pH final da solução foi de 2.8 e a produção de H,SO,, 
calculada a partir do S(6) presente na solução, foi de 58.8 kg/ton de rocha. 

Na Fig. 3 são apresentados o script do modelo de entrada (Fig. 3a) e parte do 
arquivo de saída da simulação (Fig. 3b) onde destacam-se a composição da solução final 
(Solution composition e Description of solution) e as solubilizações dos minerais do modelo 
(Phase assemblage). Na Tabela 5 é apresentada a distribuição das espécies do S(6). 
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SOLUTION 1 
temp 25 
PH 7 
pe 4 
redox pe 
units mmol/kgw 
density 1 
-water 1 fkg 
EQUILIBRIUM PHASES 1 
coz (g) -3.5 10 
Illite 00.3 
K-Feldspar 0 2.9 
02 (9) -0.7 10 
Pyrite 00.3 


SELECTED OUTPUT 1 
-file 
-reset 
-totals 
-equilibrium phases 


-gases 
END 


TITLE POTENCIAL MÁXIMO DE ACIDEZ 


teste augusto.x1s 
true 
o(o) 


K-Feldspar Illite Pyrite kaolinite 


Fe(0H)3 Gibbsite O2(g) 
Montmor-K Quartz 
oz (9) 


a 


Jarosite 


Reaction step 1. 


Usíng solution 1. 
Using pure phase assemblage 1. 


Phase SI log IAP 
coz (a) 3.50 11.33 
Illite 0.00 8.88 
K-Feldspar 0.00 0.38 
oz(g) 0.7 9 


Pyrite 


Percent error, 


Ma 
Initial 


log K(T, P) 


Tot 
Total O = 5.8 


joles in assemblage 


Final Deita 


Charge balance 
Adjusted to redox equilibrium 


Figura 3: Script de entrada da simulação (a) e da saída (b). 


Concentração Total S(6)= 6.007e-01 mol 


Distribuição das Espécies 


so,? 
AISO,); 
AISO,* 
MgSOo, 
HSO, 
FesO,* 
Fe(SO,); 
Kkso, 
Feso, 
H,so, 
KHSO, 


Molalidade 
2.61e-01 
1.29e-01 
5.36e-02 
2.08e-02 
6.52e-03 
1.02e-03 
5.33e-04 
6.89e-08 
2.08e-08 
5.82e-09 
5.74e-11 


Atividade 
3.09e-02 
8.65e-02 
3.61e-02 
2.08e-02 
4.39e-03 
6.86e-04 
3.58e-04 
4.63e-08 
2.08e-08 
5.82e-09 
5.74e-11 


Tabela 5: Distribuição das espécies na solução após ter sido alcançado o equilíbrio com a rocha na 


TESTE NAG DE ADIÇÃO ÚNICA 


O potencial de geração ácida da amostra do rejeito do BF4 foi calculado a partir 
do Teste NAG, proposto por Amira (2002). Resumidamente o teste tem como função a 
obtenção da estimativa da quantidade de ácido sulfúrico gerada em uma tonelada de 
amostra, a partir da reação de 2.5 g de amostra com peróxido de hidrogênio (15%) e sua 
titulação, até o pH 4.5, com uma solução de NaOH. Na Tabela 6 são apresentados os 


resultados obtidos. 


Descobertas e desafios nas geociências 


modelagem geoquímica. 


Capítulo 4 


65 


Amostra NAGpH o NaOH molde NAG(kg H,S0,1) 
Ai 2.59 740 0.00071 13.92 
A2 2.59 7.95 0.00080 15.58 
A3 2.61 710 0.00071 13.92 


Tabela 6: Resultados obtidos do ensaio Teste NAG. 


As amostras de rejeito do BF4 possuem um NAGpH menor que 4.5 e NAG maior que 


5. Assim, conforme a Tabela 3, os rejeitos são classificados como “Potencial de Formação 


de Ácido”. 


Os mesmos resultados foram alcançados por meio da modelagem geoquímica. 


Utilizando os 2.5 g de rocha inicialmente em contato com uma solução aquosa na qual 


foi adicionado 0.05 mol de O,, o que levou à solubilização completa da pirita e gerou uma 


solução com pH 3.89. Nessa solução foram adicionados 0.0007 mols de NaOH para que o 


pH alcançasse o valor de 4.53 (Figs. 4 e 5). Na modelagem vários minerais apresentaram 


índice de saturação superior a 1, indicando tendência à sua precipitação. Porém, nos 


experimentos não foi observada a precipitação de nenhuma fase sólida. Por essa razão na 


modelagem também não foi permitida a precipitação de minerais. 


TITLE TESTE NAG DE ADIÇÃO ÚNICA fase 1 


SOLUTION 2 
temp 25 
PH 7 
pe 4 
redox pe 
units mmol/kgw 
density 1 
-water 1 f£kg 
EQUILIBRIUM PHASES 1 
coz (g) -3.5 10 
Illite 0 0.00073 
K-Feldspar 0 0.0073 
o2(g) -0.7 0 
Pyrite 0 0.00073 
REACTION 2 
o2(g) 1 


0.05 moles in 1 steps 
SAVE solution 3 
END 
TITLE TESTE NAG DE ADIÇÃO ÚNICA - fase 2 
USE solution 3 


REACTION 3 

NaOH 1 

0.01 moles in 100 steps 
END 


Reaction 2. 


5.000e-02 moles of the following reaction have ben added: 


Relative 
Reactant moles 


o2(9) 1.00000 


Element 


Mo] 
Phase SI log IAP log K(T, P) Initial 


pH = 3.858 Charge balance 
16.703 Adjusted to redox equilibrium 


Percent error, 


Total H = 1.110507e+02 
Total O = 5.553475e+01 


b 
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Figura 4: Scripts de entrada de dados para a modelagem geoquimica fases 1 e 2 (a) e seus respectivos 
arquivos de saída (b) e (c) . Em (d) é apresentada a variação do pH em função da adição do NaOH, 
obtida pela modelagem geoquímica. 


Neutralização ácida 


B. C. RESEARCH 


O potencial de neutralização acida da amostra do rejeito do BF4 foi calculado a 
partir do teste B. C. Research (Duncan & Walden 1975). O teste consiste em uma titulação 
direta de 10 g de amostra diluída com 100 ml de água deionizada por titulação de uma 
solução padrão de H,SO, até que seja alcançado o pH 3.5. Na Tabela 7 são apresentados 
os resultados experimentais. 


MOLS DE H,SO, 


AMOSTRA VOLUME DE H,SO, (ML) ADICIONADOS PN (KG H,SO,/TON) 
Ai 3.14 0.000409268 0.4 
A2 2.42 0.000315423 0.3 
A3 3.08 0.000401447 0.4 
MÉDIA 2.88 0.000375379 0.4 
MODELAGEM 0.000361 0.4 
GEOQUIMICA 


Tabela 7: Resultado obtidos no ensaio B.C. Research e na modelagem geoquímica 


O PN apresenta uma baixa quantidade de ácido sulfúrico que pode ser neutralizada 
por tonelada da amostra de rejeito. Esse valor deve ocorrer em virtude da ausência de 
minerais capazes de neutralizar o ácido da formação na DAM, como carbonatos e 
hidróxidos, que neutralizam. 
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A baixa capacidade de neutralização é fruto do pequeno consumo de prótons, 
consequência da baixa reatividade das rochas. Em 10 g de amostra há 3*102 mol de 
microclinio, 3*10º mol de ilita e de pirita. Na modelagem geoquímica da neutralização ácida 
até o pH 3.5 por meio da adição do H,SO, (Fig. 5), replicando os experimentos, houve a 
solubilização de 1,5*10 mol de ilita e pirita e de 1,5*104 mol de microclínio (Fig. 5). Vale 
destacar que essas rochas já produzem DAM na mina Osamu Utsumi, e em virtude disso, 
sua reatividade com soluções ácidas é baixa por causa dos minerais já consumidos pela 


lixiviação prévia. 


TITLE NEUTRALIZAÇÃO ÁCIDA -  B. C. RESEARCH 
INCREMENTAL REACTIONS False 
SOLUTION 1 
temp 25 
pH 7 
pe 4 0.004 1 
redox pe 
units mmol/kgw 
density 1 os 
-water 0.1 £ kg 


0.002) 
EQUILIBRIUM PHASES 1 


mols de H2S04 adicionado 


Illite 0 1.5e-05 Dodi 
K-Feldspar 0 0.00015 
o2(g) -0.7 10 
E 0.000 ) 
Pyrite 0 1.5e-05 
Quartz 0 Es 
REACTION 1 
H2504 0.1 
0.0103 moles in 100 steps b 
END 
a 
e e e e e e mm Phase assemblage----——————————————————————————— —— 
Moles in assemblage 
Phase SI log IAP log K(T, P) Initial Final Delta 
Illite =5.53 3.35 8.ss 1.500e-05 oO -1.500e-05 
K-Feldspar -3.88 -4.26 -0.38 1.500e-04 o -1.500e-04 
o2(g) -0.70 =3.59 -2.89 1.000e+01 1.000e+01 -8.175e-05 
Pyrite -226.94 -251.64 -24.70 1.500e-05 0 -1.500e-05 
Quartz 0.00 -4.03 -4.03 0.000e+00 4.931e-04 4.931e-04 


Cc 


Figura 5: Script de entrada do modelo geoquímico (a) e a solubilização dos minerais para alcançar o pH 
3.5 (c). Em (b) é apresentada a alteração do pH da solução em função da adição do H2S04 0,0108 M. 


CONCLUSÕES 


Os métodos de avaliação de geração e neutralização ácida permitem uma simples 
avaliação sobre a drenagem ácida de mina. O potencial máximo de acidez (PMA) obteve 
um valor 55.51 kg de H,SO,/ton de amostra e se mostrou eficiente na quantificação do total 
de ácido gerado na DAM. O teste NAG quantificou 14.47 kg de H,SO (ton de amostra. Esse 
valor é menor que o PMA devido ao tempo de ensaio não ser suficiente para toda amostra 
reagir. 
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A falta de minerais capazes de consumirem prótons em suas solubilizações, 
conforme resultados dos testes do método B.C. Research, comprova o baixo potencial de 
neutralização destas rochas do BF4 da mina de Osamu Utsumi. 

Por meio da modelagem geoquímica foi possível reproduzir os ensaios laboratoriais 
e predizer as capacidades de geração e neutralização da acidez. Além disto, por meio da 
modelagem foi possível estimar a quantidade solubilizada de minerais, responsável pela 
baixa capacidade de neutralização das rochas. 

Desta forma, a partir da determinação mineralógica com quantificação por modelagem 
Rietveld, é possível se estimar as capacidades de geração e de neutralização de DAM 


antes mesmo de se realizar os experimentos utilizando-se a modelagem geoquímica. 
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